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SECTION  1 

INTRODUCTION 


The  Aeronutronlc  program  on  the  detonation  behavior  of  propel* 
lanta  is  directed  to  the  development  of  mathematical  techniques  by  which 
such  behavior  can  be  determined  solely  from  basic  hydrodynamic  princi¬ 
ples  and  properties  of  the  material.  Particular  intereat  ia  in  the  pre¬ 
diction  of:  (1)  initiation  behavior,  and  (2)  energy  yield  and  the  sta¬ 
bility  of  detonation,  as  a  function  of  type  of  initiation,  material 
parameters  and  geometry. 

The  approach  which  has  been  adopted  Involves  representation 
of  the  detonation  process  in  terms  of  a  set  of  basic  equations  consisting 
of  the  Navier-Stokes  conservation  equations,  the  appropriate  equations  of 
state,  and  the  rate  equations  which  describe  the  release  of  chemical 
energy.  This  system  of  equations  is  integrated  numerically  in  time  and 
space  making  use  of  electronic  computer  techniques  to  obtain  the  history 
of  any  defined  initial  disturbance.  Computation  of  such  histories  for 
various  configurations  and  initial  conditions  will  define  detonation 
behavior. 


The  prosecution  of  the  program  has  been  in  several  phases 
worked  on  more  or  less  concurrently:  (1)  development  of  the  proper  form 
of  the  basic  equations,  (2)  development  of  suitable  mathematical  techniques 
for  solution  of  the  basic  equations,  (3)  development  of  the  necessary 
input  data  on  pertinent  propellant  materials,  (4)  performance  of  calcula¬ 
tions  on  explosive  charges  of  experimentally  known  behavior,  and  compari¬ 
son  of  computed  and  experimental  results,  and  (5)  calculation  of  the 
behavior  of  propellant  grain  configurations  of  particular  interest  to 
the  propulsion  field  to  determine  dependence  of  behavior  upon  grain 
design  parameters. 
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Phases  1  and  2  of  this  program  have  been  completed,  and  some 
work  has  been  carried  out  on  phases  3  and  4.  The  basic  equations  have 
been  written  in  a  form  which  now  appears  to  correctly  incorporate  at 
least  those  parameters  and  material  characteristics  which  are  of  great* 
est  significance  to  the  detonation  process.  Mathematical  techniques  and 
associated  computer  programs  have  been  written  and  checked  out  for  solving 
the  basic  equations  in  one  or  two  space  dimensions  and  for  pressure  init¬ 
iation.  The  problem  of  thermally  initiated  detonation  has  been  extensively 
studied  with  the  conclusion  that  it  is  beyond  the  capability  of  the  present 
generation  of  computers.  Experimental  data  have  been  obtained  on  the 
equation  of  state  parameters  for  the  unreacted  solid  phase.  Similar  data 
for  the  reacted  gaseous  products  have  been  obtained  from  the  detonation 
literature. 


As  of  the  present  time,  considerable  success  has  been  achieved 
in  demonstrating,  on  the  part  of  theoretically  computed  detonations, 
general  agreement  with  experimentally  observed  behavior.  This  has  served 
to  verify  the  validity  of  the  theoretical  methods,  and  in  doing  so,  it 
has  contributed  substantially  to  our  understanding  of  the  detonation 
process  since  the  factors  responsible  for  the  various  aspects  of  detona¬ 
tion  behavior  are  readily  seen  from  the  basic  equations.  For  example, 
the  initiation  of  a  detonation  wave  is  found  to  be  dependent  upon  the 
equation  of  state  of  the  unreacted  material,  the  kinetics  of  the  chemical 
reactions  and  the  heat  of  reaction.  The  concept  Of  the  ignition  zone 
which  developed  from  this  work  and  its  relation  to  the  hot  spot  theory 
of  Initiation  as  advanced  by  Bowden  and  Yoffe  have  made  clear  the  details 
of  such  dependence.  Wave  growth  following  initiation  has  been  shown  to 
be  in  accord  with  the  classical  grain  burning  model  of  Eyrlng  and 
coworkers;  however,  in  some  respects  a  point  ignition  internal  burning 
model  appears  to  be  more  realistic.  Finally,  the  phenomena  of  low  order 
detonation  which  heretofore  was  not  capable  of  clear  theoretical  explana¬ 
tion,  has  been  shown  to  be  a  normal  consequence  of  certain  combinations 
of  the  above  processes. 

In  addition  to  its  success  in  demonstrating  and  explaining  the 
general  aspects  of  detonaticxi  behavior,  the  theoretical  method  has  also 
been  used  to  ccmpute  the  quantitative  initiation  behavior  of  real  mater¬ 
ials.  Such  efforts  are  severely  handicapped  by  a  lack  of  the  material 
data  needed  for  input  to  the  problem.  However,  in  the  case  of  two  pro¬ 
pellants  for  which  such  data  were  available,  calculation  of  minimum 
pressures  for  initiation  shewed  excellent  agreement  with  experiment. 

Remaining  work  to  be  done  involves  (1)  the  procurement  of  needed 
additional  data  on  the  various  equations  of  state  and  on  the  parameters 
appearing  in  the  equations  which  describe  the  rate  of  release  of  chemical 
energy,  (2)  carrying  out  additional  calculations  on  the  detonation  behavior 
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of  propellants  and  high  explosives  which  can  be  compared  with  experimental 
observations  for  the  purpose  of  further  definition  of  the  validity  of  the 
mathematical  methods,  and  (3)  performance  of  calculations  which  are  of 
application  to  the  development  and  practical  evaluation  of  new  propellant 
systems. 
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SECTION  2 

THE  PHYSICAL  MODEL 


2.1  General  Considerations 

A  detonation  is  represented  as  a  hydrodynamic  shock  in  which 
chemical  reaction  takes  place  in  the  expansion  zone  behind  the  shock 
front.  The  expansion  converts  part  of  the  heat  of  reaction  into  work 
and  the  combination  of  this  added  work  with  that  lost  in  the  shock  by 
dissipation  effects  corresponding  to  the  entropy  increase,  radiation, 
thermal  conduction,  lateral  expansion,  etc.  is  the  condition  which 
determines  the  wave  history.  When  the  work  loss  exceeds  the  work  added 
by  reaction  the  wave  fades;  when  the  work  added  by  reaction  exceeds  the 
loss,  the  wave  grows,  when  the  two  are  the  same,  the  wave  is  at  steady 
state. 


A  process  such  as  that  described  can  be  represented  mathe- 
matically  by  a  set  of  basic  equations  consisting  of  (1)  the  equations 
of  conservation  of  mass,  momentum  and  energy,  (2)  the  equations  describ¬ 
ing  the  rate  of  release  of  chemical  energy,  and  (3)  the  equations  of 
state  of  the  various  material  phases  taking  part  in  the  process.  Three 
requirements  exist:  first,  the  equations  must  be  stated  in  a  form 
capable  of  representing  the  model  in  a  fashion  sufficiently  detailed  for 
the  requirements  of  the  treatment;  second,  the  parameters  which  appear 
in  the  equations  representing  the  chemical  and  physical  properties  of 
this  material  must  be  known,  and  third,  techniques  for  obtaining  solu¬ 
tions  to  the  system  of  expiration  must  be  available. 


2.2  The  Conservation  Equations 

The  numerical  techniques  (to  be  discussed  later)  used  in  the 
treatment  of  detonation  problems  are  most  conveniently  applied  when  the 
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conservation  equations  are  expressed  in  Lagrange  form.  If  the  medium  is 
a  fluid,  only  pressure  and  viscous  forces  need  be  considered,  otherwise 
additional  forces  due  to  mechanical  rigidity  are  present.  Since  pressures 
of  interest  to  detonation  phenomena  generally  exceed  the  mechanical 
strength  of  materials  by  a  wide  margin,  the  assumption  will  be  made  that 
all  materials  act  as  fluids.  Viscous  forces  affect  detonation  behavior 
only  in  terms  of  the  work  dissipation  effects  occurring  at  the  shock 
front.  It  is  convenient  to  include  such  forces  in  an  artificial  viscosity 
term,  "q",  original  devised  by  von  Newman  and  Richtmeyer*  which  is 
designed  not  only  to  provide  the  correct  work  dissipation  in  a  shock,  but 
also  to  convert  the  mathematical  form  of  the  shock  from  a  discontinuity, 
which  it  assumes  in  a  nonviscous  medium  to  a  continuous  function  in  which 
rates  of  change  are  small  enough  to  permit  numerical  integration.  This 
greatly  simplified  treatment  of  the  problem  in  the  shock  region. 


With  these  considerations,  the  conservation  expirations  are 
the  following: 

Mass  |-TT"  +  *M  ,v|p  *  *  *P  (1) 


Momentum  pi  -^7  +  p  •  v)  p  -  -  V  (  p+q) 
Energy  juvje--|p+q  V-  p 


where  q  ■ 


aT+  * 


P  +  q  v.  p 


if  d  u/dx<  0 
if  d  u/dx^  0 


and  X  are  an  arbitrary  constants  which  define  the  width  of  the 


region  of  space  over  which  the  shock  zone  is  spread  and  the  rate  of  change 
of  the  variables  in  this  region. 


2.3  Equation  of  State 

Reaction  in  the  detonation  of  a  condensed  explosive  involves 
the  transformation  solid  (or  liquid)  — r  gas.  As  previously  mentioned, 
it  is  assumed  that  under  detonation  conditions,  both  solids  and  liquids 
behave  a  condensed  fluids,  so  that  only  two  phases  need  to  be  considered 
which  are  empirically  designated  as  solid  and  gas.  The  detonation  process 
therefore  is  postulated  to  initiate  in  a  solid  phase,  and  terminate  in  a 
gaseous  phase.  The  situation  in  between  is  dependent  upon  the  nature  of 
the  reaction  process.  This  is  discussed  in  the  following  section  where 
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it  will  be  shown  that  the  most  reasonable  concept  of  the  reaction  process, 
due  to  Eyring  and  co-workers^  is  consistent  with  the  concept  of  the 
partially  reacted  material  as  a  physical  mixture  of  the  solid  and  gaseous 
phases . 


If  these  phases  are  physically  different  the  difference  will 
appear  as  a  difference  in  their  respective  equations  of  state.  In  accord¬ 
ance  with  the  conception  of  the  solid  phase  as  a  fluid,  any  such  distinction 
should  become  small  and  the  two  equations  of  state  should  tend  to  become 
identical  at  high  pressure.  At  low  pressure,  on  the  other  hand,  they 
should  reduce  to  the  conventional  gas  law,  and  a  solid  equation  of  state 
which  properly  describes  the  important  aspects  of  the  behavior  of  solids 
at  ambient  conditions.  It  is  desirable  to  consider  the  form  of  equations 
which  possess  these  characteristics. 

a-  Equation  of  State-Gas  Phase 

The  equation  of  state  of  the  gaseous  products  of  a  detonation 
process  appears  in  the  steady  state  theory  of  detonation,  and  comparison 
of  the  ideal  detonation  velocity  predicted  by  steady  state  theory  with 
experimental  value  provides  some  insight  into  the  nature  of  this  equation. 
The  general  form  seems  to  be  that  of  the  perfect  gas  law  with  one  arbi¬ 
trary  variable  parameter.  Specifically  the  co-volume  equation  of  state 
with  co-volume  term  as  a  function  of  volume:^ 


P  V  -  a(v)] 

|  -  RT 

(5) 

the  form  with  a  pressure  dependent  parameter^ 

pv  «  nRT  + 

a(p) 

(6) 

the  Kistiakowsky-Wilson  equation^ 

„„  ,  nRT 

i 

pv 

1  +  xe  1 

Kf) 

with  x  -  £(1^  kj/T^vM 

(8) 

and  the  Boltzman  equation  of  stated  (in  various  modifications) 


dv  b«  b03 

-  i  .  +  0.625  —2 1  +  0.2869  — §-  (9) 

nRT  v  v  v 

have  all  been  used  with  some  success  in  computing  various  aspects  of 
steady  state  behavior  which  agree  well  with  experiment.  In  none  of  these 
cases  has  it  been  possible  to  determine  the  variable  parameter  from  inde¬ 
pendent  data,  a  fact  which  makes  difficult  any  Independent  choice  between 
them. 
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Some  work  which  is  believed  to  have  an  Important  bearing  on 
this  problem  concerns  studies  on  the  equation  of  state  of  solids  and 
liquids  at  detonation  pressures. ^ •'  It  suggests  that  they  are  well  rep¬ 
resented  by  a  simple  perfect  gas  equation  of  state  in  the  form: 


pv  -  (7-1)  e 


(10) 


if  a  value  y  -3  is  used.  In  accord  with  the  principle,  cited  above,  that 
solids  and  gases  should  possess  similar  equations  of  state  at  high 
pressure,  it  is  deduced  that  the  perfect  gas  equation  of  state  is  also 
valid  for  gases  at  high  pressure  provided  an  appropriate  value  of  is 
used.  These  concepts  suggest  an  equation  of  state  of  the  form  of  (10) 
for  the  gas  phase  in  which  7  is  the  variable  parameter  having  values  in 
the  range  of  1.3  for  perfect  gases  at  1  atmosphere  to  7~3  at  very  high 
pressures. 


Such  an  equation  of  state  has  been  examined  by  calculating  the 
range  of  values  needed  to  be  assumed  by  7  in  order  for  equation  (1),  when 
used  in  conjunction  with  steady  state  theory,  to  provide  detonation 
velocities  in  accord  with  experiment.  Writing  5  in  a  form  analagous  to 
(10),  i.e. 


g 


Pv„ 


p  v  -a(v  )  -  (yo-l)e 


g' 


7„  -  1  + 


g 


£2 


one  obtains: 


i\+& 

e 

g 


7 


g 


(11) 


(12) 


3 

By  means  of  studies  on  a  series  of  explosives,  M.A.  Cook  exper¬ 
imentally  determined  the  function  a  (v).  By  comparison  of  (5)  with  the 
virial  equation  of  state,  it  is  found  that  a  can  be  satisfactorily  repre¬ 
sented  by  the  expression 


2 

<*  «  v(l  -  -s - - - j”)  , 

v  +  bv  +  b^ 


(13) 


in  which  the  first  three  terms  of  the  virial  expansion  are  used. 
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Combination  of  (12)  and  (13)  leads  to: 

7g  "  70+-£“  <l  +  £->  <70  -  D-  04) 

g  g 


The  validity  of  expression  (14)  is  shown  in  Table  I,  where 
experimental  values  of  a  as  obtained  by  Cook  are  to  be  compared  with 
values  computed  from  (14)  at  a  series  of  values  of  v^.  The  agreement 

shown  is  surprisingly  good,  in  view  of  the  fact  that  application  of  (14) 

is  extended  to  the  condition  where  ->1  and  the  virlal  equation  of  state 
no  longer  converges.  v 


Some  interpretation  as  to  the  meaning  of  the  parameter  7  is 
desirable.  As  it  is  used  in  equation  5,  it  can  be  shown  to  be  defined  by: 


(15) 


For  a  perfect  gas  ( o)  this  reduces  to: 


7  " 


din  p  | 

d  lnP  Is 


(16) 


which  can  be  shown  to  be  equivalent  to  the  familiar  relation: 

C 


(17) 


Any  departure  of  the  observed  y  from  (17)  can  be  interpreted  as  the 
appearance  of  a  significant  contribution  by  the  term 


1/p 


dpa 

dv  S 


Analagous  lines  of  reasoning  can  be  followed  for  the  other  equations  of 
state,  6,  7,  and  9. 

b.  Equation  of  State  -  Solid  Phase 

The  body  of  work  which  has  become  available  on  the  equation  of 
state  of  condensed  materials  in  recent  years  is  much  too  extensive  to  be 
discussed  comprehensively  here.  The  most  pertinent  information  appears 
to  derive  from  studies  on  shock  propagation? «?and  these  lead  to  observa¬ 
tions  that  at  pressure  in  the  detonation  range,  a  simple,  perfect  gas 
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CALCULATION  OF  a  FROM  EQUATION  (13)  and  7  FROM  EQUATION  (14) 
COMPARISON  WITH  EXPERIMENTAL  DATA 


v(cm^) 

a(exp) 

a(calc.)* 

I? 

1.25 

0.75 

0.751 

1.75 

1.0 

0.68 

0.676 

1.93 

0.837 

0.62 

0.611 

2.12 

0.714 

0.56 

0.555 

2.36 

0.625 

0.51 

0.507 

2.61 

0.555 

0.47 

0.466 

2.89 

0.50 

0.43 

0.431 

3. 19 

0.455 

0.39 

0.400 

3.52 

0.19 

0.17 

0.185 

11.7 

*  b  -  1.03,  7  -  1.3 
o 
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equation  of  state 


pv  -  (7  -l)e 

where  7  has  a  value  of  about  3,  la  applicable  within  experimental 
uncertainty.  At  lower  pressures,  this  expression  is  clearly  not  appli¬ 
cable  since  it  does  not  provide  for  a  finite  volume  at  zero  pressure. 
This  suggests  the  addition  of  additional  parameters  P  and  *  to  give 
an  equation  of  the  following  form: 


(p  +  £  )  v  ■  (  7S  -  l)e  -  e  (18) 

Equation  (18)  is  of  the  form  of  an  equation  of  state  deduced 
for  solids  by  Gruneisen,  based  upon  the  Einstein  concept  of  the  solid 
phase  as  a  collection  of  harmonic  oscillators.  If  the  possible  frequen¬ 
cies  are  v  ,  the  equation  obtained  from  such  a  model  is  the  following: 


(19) 


where 


Here  the 
and  0  is 
given  by 


summed  terms  represent  the  internal  thermal  energy  of  the  material 
a  potential  function  such  that  the  Internal  potential  energy  is 


(21) 


For  present  purposes,  Internal  thermal  energy  is  meant  to  include  the  kinetic 
and  potential  energy  associated  with  atomic  and  molecular  vibration,  and 
internal  potential  energy  represents  that  energy  due  to  atomic  and  molecular 
distortion. 


At  high  temperature  the  summation  in  equation 

«  T?  i  • 5 

where  £  is  the  zero  point  energy.  If  all  of  the  'y  are 
of  (22)  in  (19)  gives:  a 


(19)  reduces  to 
(22) 

equal,  substitution 
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(p+^-)V-'VG  (3NkT’S)  <23> 

The  oscillator  model  upon  which  equation  (23)  is  based  is  closely 
approximated  only  by  a  perfect  single  crystal.  However,  the  form  of  the 
equation  is  a  consequence  of  the  condition  that  the  potential  and  thermal 
forms  of  internal  energy  are  respectively  dependent  upon  the  volume  only 
and  upon  the  temperature  only,  and  aretherefore  separable  into  different 
terms.  If  it  is  assumed  that  however  poorly  the  simple  oscillator  model 
may  represent  reality  in  other  respects,  this  separation  of  thermal  and 
potential  energy  is  valid  for  real  materials,  then  the  ..equation  is  of 
general  application,  provided  only  that  the  functions  ^  and  7  ,  which 
now  must  be  regarded  as  somewhat  empirical  in  nature,  are  known. 


In  treating  detonation  behavior,  it  is  more  convenient  to 
express  the  equation  of  state  in  terms  of  total  energy,  since  the  conser¬ 
vation  equations  refer  to  total  energy.  Adding  equation  (21)  to  both 
sides  of  (23),  one  has  ultimately 


‘^Ge+7G 


[3  nkTo+|] 


-  —'0  "  L  ■* 

This  is  now  the  form  of  the  solid  equation  of  state  presently  in 


(24) 

use  where 


(25) 

(26) 


The  quantities  j3 1  and  £'  can  be  expressed  in  terms  of  their 
values  (3  and  £  at  ambient  conditions  as  follows: 


?  "  ^(U’V  v) 

£ '  -  £  +  i  (7g) 

£ 

where  fami  *  o 

£  e 

at  ambient  conditions. 


(27) 

(28) 
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One  obtains 


‘V1’ 

e  -(f§/3v§  +  f86)/(78"l) 

jk\f  v  +v 

1  I  s  s 

(34) 


An  equation  for  pressure  during  the  ignition  phase  is  obtained 
through  use  of  expression  (18) •  Upon  differentiation  we  obtain 

(p  +  (S  )  dv  +  vdp  ■  (7  -1)  de  +  ed7 

combination  with 

de  ■  - (p  +q)  dv 

and 


illhf 

v  I  5x1 


+  C.v 


Du 


gives,  finally,  to  a  good  approximation: 


(P 


[P  +  £  -  e/v 

P+£/7  + 


constant  (34a) 


2.4  The  Chemical  Rate  Equations 
a.  General  Model 


Reaction  in  a  detonation  wave  is  divided  into  two  distinct 
phases;  an  ignition  phase  in  which  follows  an  Arrhenius  rate  law,  and  a 
subsequent  burning  phase  based  upon  the  Eyring  grain  burning  model. 

These  two  phases  occur  in  succession,  and  involve  separate  portions  of 
the  charge  material.  The  fraction  of  the  charge  Included  in  each  is  one 
of  the  parameters  of  the  material  and  must  be  supplied  as  input  to  a 
given  problem. 

b.  Ignition  Phase 

The  portion  of  the  charge  which  takes  part  in  the  ignition 
phase  is  based  upon  the  concept ,  somewhat  analagous  to  the  hot  spot 
theory  of  initiation,  that  only  certain  discrete  portions  of  the  charge 
are  caused  to  react  by  the  event  which  initiates  detonation.  This  is 
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Substituting  and  transposing 


(P  +  j8)  V  -7  Ge  +  e  +  /£  -v  / 

(29) 

Setting 

(30) 

One  has  finally  the  previously  deduced  expression 

(p  +  A)  v  +  (y-l)e  +  e 

(18) 

When  £  and  e  are  constants  obtainable  from  the  properties  of  the  material 
at  ambient  conditions  and  CY-1)  la  a  complicated  function  of  7  ,  v 

av  G 

and  e,  equal  to 7_  at  ambient  conditions.  Equation  18  was  used  as  the 
G 

solid  equation  of  state  throughout  the  present  work,  with  (yg  -1)  taken 
as  an  experimentally  determined  variable  parameter. 


Although  only  fragmentary  data  are  available  from  which  to 
determine  the  value  of'V'-l  for  solids  at  high  pressure,  it  is  a  logical 
deduction  that  gaseous  and  condensed  phases  are  essentially  indisting¬ 
uishable  at  detonation  pressures.  Accordingly  the  assumption  was  made 
that  the  information  discussed  in  the  previous  section  on  the  value  of 
7g  for  the  reaction  products  of  detonation  could  be  applied  directly  to 

define  7  for  the  solid  phase.  Accordingly ,7g  was  computed  by  an 

expression  analagous  to  14,  namely 


7 

8 


y  + 

0 


b_ 

v 

8 


+ 


b_ 

vg 


7 

o 


(31) 


Equations  (18)  and  (31)  were  the  ones  used  throughout  the  program  to 
represent  the  solid  phase. 

In  regions  where  the  charge  consists  of  a  half  reacted  mixture 
of  solid  reactant  and  gaseous  products,  it  is  necessary  to  have  an 
equation  giving  the  pressure  in  terms  of  the  relative  proportions  of 
solid  and  gas,  and  the  energy  density.  This  is  found  by  eliminting  eg , 

e  and  v  between  equations  10  and  18  and 
8  8 

*  -7,  v#  +  (1  -  7,)  vg 
*-7,e#  +  (l-7,)eg 
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presumed  to  be  a  consequence  of  an  effect,  either  of  such  areas  being 
raised  by  the  initiation  event  to  a  higher  temperature  than  their  sur¬ 
roundings  (hot  spot  theory)  or  of  being  more  reactive.  Reaction  in  these 
regions  is  assumed  to  be  represented  by  an  Arrhenim  function  of  the  form 

St  “  (1*f)  exp(ir)  (35) 

As  discussed  in  the  previous  section  on  equation  of  state,  expression  is 
in  terms  of  total  energy  density,  e,  rather  than  temperature.  If  reaction 
is  dependent  only  upon  thermal  energy  density,  then  the  form  of  35  in¬ 
volves  the  tacit  assumption  that  all  energy  changes  are  thermal.  This 
is  likely  to  be  erroneous  to  some  degree  at  detonation  pressures. 

TWo  physical  conceptions  of  the  ignition  region  of  the  charge 
are  possible  which  have  significant  implications  in  terms  of  the  follow¬ 
ing  reaction  process  to  be  discussed  in  the  following  section.  The  first 
is  based  upon  a  grain  structure  for  the  detonable  material  and  conceives 
the  ignition  region  to  be  a  thin  skin  of  material  on  the  grain  surface. 
Material  in  this  skin  is  assumed  to  be  more  reactive  than  that  in  the 
middle  of  the  grain,  and  to  react  uniformly  throughout  its  volume  the 
latter  being  required  in  order  for  it  to  exhibit  a  reaction  rate  accord¬ 
ing  to  35.  This  assumption  can  be  justified  on  the  basis  of  arguments 
that  reaction  at  the  interior  of  a  lattice  structure  involves  creation 
of  a  discontinuity,  a  highly  endothermic  process.  Furthermore,  reactive 
solids  are  normally  observed  to  react  at  interfaces  or  discontinuities. 

The  second  concept  requires  no  formal  grain  structure,  but 
visualizes  the  ignition  region  to  consist  of  a  large  number  of  more  or 
less  regularly  spaced  regions  of  finite  volume  located  throughout  an 
otherwise  continuous  material.  These  regions  are  postulated  to  be  diff¬ 
erent  from  contiguous  material  such  that  when  an  initating  event  occurs, 
the  material  in  these  regions  reacts  first  and  according  to  the  rate  law 
(35). 

9  10 

They  correspond  to  the  hot  spot  concept  of  Bowden  and  Yoffee,  ’ 
explained  by  those  authors  as  due  to  such  mechanisms  as  contact  between 
grit  particles,  compression  of  small  voids  and  the  like.  It  is  of 
interest  to  consider  one  of  these  in  some  detail,  namely,  Ignition 
regions  due  to  the  pressure  of  small  voids  or  gas  bubbles. 

A  number  of  Investigators  have  considered  the  possibility  that 
the  sensitizing  effect  of  a  porous  structure  upon  explosives  was  due  to 
the  presence  of  gas  in  the  pores,  which,  when  adiabatically  compressed 
by  a  transient  pressure,  would  be  heated  to  temperature  above  the  ignition 
point  of  the  surrounding  material.  There  is  no  question  but  that  occluded 
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presumed  to  be  a  consequence  of  an  effect,  either  of  such  areas  being 
raised  by  the  initiation  event  to  a  higher  temperature  than  their  sur¬ 
roundings  (hot  spot  theory)  or  of  being  more  reactive.  Reaction  in  these 
regions  is  assumed  to  be  represented  by  an  Arrhenim  function  of  the  form 

f  ■  <‘-f>  «p(ir]  <35> 

As  discussed  in  the  previous  section  on  equation  of  state,  expression  is 
in  terms  of  total  energy  density,  e,  rather  than  temperature.  If  reaction 
is  dependent  only  upon  thermal  energy  density,  then  the  form  of  35  in¬ 
volves  the  tacit  assumption  that  all  energy  changes  are  thermal.  This 
is  likely  to  be  erroneous  to  some  degree  at  detonation  pressures. 

Two  physical  conceptions  of  the  ignition  region  of  the  charge 
are  possible  which  have  significant  implications  in  terms  of  the  follow¬ 
ing  reaction  process  to  be  discussed  in  the  following  section.  The  first 
is  based  upon  a  grain  structure  for  the  detonable  material  and  conceives 
the  ignition  region  to  be  a  thin  skin  of  material  on  the  grain  surface. 
Material  in  this  skin  is  assumed  to  be  more  reactive  than  that  in  the 
middle  of  the  grain,  and  to  react  uniformly  throughout  its  volume  the 
latter  being  required  in  order  for  it  to  exhibit  a  reaction  rate  accord¬ 
ing  to  35.  This  assumption  can  be  justified  on  the  basis  of  arguments 
that  reaction  at  the  interior  of  a  lattice  structure  involves  creation 
of  a  discontinuity,  a  highly  endothermic  process.  Furthermore,  reactive 
solids  are  normally  observed  to  react  at  interfaces  or  discontinuities. 

The  second  concept  requires  no  formal  grain  structure,  but 
visualizes  the  ignition  region  to  consist  of  a  large  number  of  more  or 
less  regularly  spaced  regions  of  finite  volume  located  throughout  an 
otherwise  continuous  material.  These  regions  are  postulated  to  be  diff¬ 
erent  from  contiguous  material  such  that  when  an  initating  event  occurs, 
the  material  in  these  regions  reacts  first  and  according  to  the  rate  law 
(35). 

9  10 

They  correspond  to  the  hot  spot  concept  of  Bowden  and  Yoffee,  ’ 
explained  by  those  authors  as  due  to  such  mechanisms  as  contact  between 
grit  particles,  compression  of  small  voids  and  the  like.  It  is  of 
interest  to  consider  one  of  these  in  some  detail,  namely,  Ignition 
regions  due  to  the  pressure  of  small  voids  or  gas  bubbles. 

A  number  of  investigators  have  considered  the  possibility  that 
the  sensitizing  effect  of  a  porous  structure  upon  explosives  was  due  to 
the  presence  of  gas  in  the  pores,  which,  when  adiabatically  compressed 
by  a  transient  pressure,  would  be  heated  to  temperature  above  the  ignition 
point  of  the  surrounding  material.  There  is  no  question  but  that  occluded 
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gas  could  be  heated  to  a  very  high  temperature  by  such  a  mechanism.  How¬ 
ever,  analysis  of  this  system  by  Johansson^  has  raised  some  doubts  as 
to  whether  the  quantity  of  heat  generated  in  this  way  would  be  sufficient 
to  cause  ignition.  He  found  it  necessary  to  postulate  the  presence  in 
the  pores  of  a  vapor  capable  of  exothermic  decomposition  so  that  addition¬ 
al  heat  could  be  accounted  for. 

An  alternative  mechanism  exists,  however,  by  which  additional 
heat  can  be  derived  from  pore  collapse.  Two  factors  are  of  importance 
in  defining  this  mechanism. 

First,  increased  porosity  is  accompanied  by  an  increase  in  the 
effective  specific  volume  and  compressibility.  This  results  in  an  in¬ 
crease  in  the  energy  deposited  in  the  material  by  a  shock  or  pressure 
transient.  Second,  this  extra  energy  is  not  deposited  uniformly  through¬ 
out  the  material.  It  appears  primarily  as  a  consequence  of  high-speed 
viscous  flow  and  other  dissipative  effects  which  accompany  collapse  of 
the  voids  as  well  as  of  the  adiabatic  or  shock  heating  of  occluded  gases 
retained  by  the  voids.  The  excess  energy  is  therefore  concentrated  in 
regions  of  the  material  contiguous  to  the  voids,  which  become  the  so- 
called  hot  spots. 

The  present  mathematical  methods  are  not  capable  of  representing 
the  charge  material  with  sufficient  resolution  to  permit  treatment  of  the 
formation  of  such  hot  spots.  However,  the  effect  can  be  simulated  by 
adjusting  the  value  of  the  activation  energy  used  in  the  ignition  zone. 


It  can  be  shown  that  for  material  which  obeys  an  equation  of 
state  18,  the  energy  deposited  by  a  heavy  shock  (p^pQ)  is  the  following: 


e 

s 


1/2  pv 


pVvim_pM 
P(7S  +  1)  +  2  ^  J 


(36) 


Neglecting^  as  small  compared  to  e,  one  has 
e  =  kv  , 


where 


k  =  1/2  p 


P 


P(7g  +  l)  +  2  £ 


(37) 
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The  increase  in  energy  Ae  deposited  for  a  given  increase  in  initial 
volume.Av,  then  becomes 

A  e  ■  kAv 

s 

If  all  of  this  increase  is  deposited  in  a  region  containing  a  fraction, 
Fp  of  the  material,  the  increase  in  energy  of  this  region  becomes 

A  es  “  [k/Fl]AV  (38) 

If  one  now  assumes  that  energy  generated  because  of  the  collapse  of  a 
void  is  deposited  in  a  volume  of  material  approximately  equal  to  the 
void  volume,  then  one  has 

Ae  =  (k/F,)F,  v  *  kv  .  (39) 

s'  V  loo  '  ' 

Thus  the  energy  density  in  a  hot  spot  created  by  the  collapse  of  a  void 
will  always  be  approximately  twice  the  energy  of  the  bulk  material. 
Simulation  of  a  "hot  spot"  ignition  region  arising  from  voids  in  the 
material  then  becomes  possible  by  use  of  an  activation  energy  one-half 
that  of  the  grain-burning  reaction  for  the  ignition  process.  It  also 
follows  that  the  term  F^,  previously  used  to  represent  the  fraction  of 
material  in  the  ignition  zone,  becomes  equal  to  the  fraction  of  the 
specific  volume  which  consists  of  voids. 


c.  Burning  Phase  -  Homogeneous  Materials 


When  material  in  the  ignition  region  has  reacted  completely, 
further  reaction  is  postulated  to  take  place  via  a  surface  regressive 
process  according  to  the  Eyring  grain-burning  concept.  It  was  shown  by 
Eyring^  that  this  led  to  a  rate  expression  of  the  form 


df 

dt 


2/3 

(1-f) 


C  A 
v 

Re 

z 


(40) 


where  ez/C  was  taken  to  be  the  explosion  temperature.  For  the  present 
application  two  modifications  are  made  in  the  equation. 


The  first  is  based  upon  the  argument  that  at  low  pressures, 
equation  40  should  approach  conventional  expressions  for  the  burning  rate 
of  a  solid  combustible  material  which  means  that  it  should  become  propor¬ 
tional  to  some  power  of  the  pressure.  Since  potential  applications  of 
the  present  treatment  include  low  pressure  processes,  a  term  was  added 
to  make  40  conform  to  such  behavior. 
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The  second  modification  is  based  upon  the  fact  that  energy 
density  at  the  regressing  grain  surface  is,  in  reality,  equal  to  the 
energy  in  the  solid  prior  to  reaction,  plus  the  heat  of  reaction,  less 
any  work  done  in  expansion  during  reaction.  It  is  assumed  that  the  sur¬ 
face  regressive  reaction  is  analagous  to  a  deflagration,  i.e. ,  a  constant 
pressure  process.  It  is  therefore  obtained  that 


With  these  two  changes,  the  form  of  grain-burning  equation 
used  in  this  work  is  the  following 


df 

dt 


(1-f) 


2/3 


aB  j> _ 

Z  p/W  +  1 


exp 


C  A 
v 

R(es 


(41) 


In  cases  where  the  ignition  region  is  taken  to  be  a  point  rather  than  a 
grain  surface  an  analagous  expression  is  obtained  by  setting 

(1-f)  *  f 


and  defining  Z  as  one  half  the  distance  between  centers, 
equation  becomes 


df 

dt 


2/3 


aB 


ZP  + 

w 


exp 


7J 


The  corresponding 


(42) 


Integration  of  equations  41  and  42  requires  knowledge  of  es  and 
7  ,  terms  which  are  variable  throughout  the  problem.  *y  is  obtained  as 
a^function  of  volume  from  equation  14.  Development  of  i  suitable  proce¬ 
dure  for  determining  e  represents  a  problem  of  significantly  greater 
complexity. 

It  was  shown  by  Eyring  and  co-workers  that  the  high  regression 
rates  of  grain  surfaces  during  the  grain-burning  process  effectively 
prevents  conduction  of  heat  into  the  unburned  solid.  Accordingly  only 
two  processes  exist  by  which  the  energy  of  the  solid  can  change;  first 
by  overall  compressions  or  expansions  due  to  hydrodynamic  motion,  and 
second,  by  an  effect  analagous  to  the  burning  of  a  solid  material  in  a 
completely  filled  bomb.  The  latter  is  shown  by  the  diagram  of  figure  1. 
The  top  diagram  (a)  shows  the  filled  bomb  before  any  material  has  burned. 
The  cross-hatched  area  at  the  left  represents  material  which  is  about  to 
react,  i.e.  the  reaction  zone.  The  second  diagram  (b)  represents  the 
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situation  after  this  material  has  reacted.  The  reacted  material  now 
occupies  a  greater  volume  because  it  obeys  a  different  equation  of  state 
and  because  of  the  released  energy  of  reaction  which  it  contains.  The 
remaining  solid  has  been  compressed  and  is  at  a  higher  energy  density 
because  of  the  compressive  work  performed  upon  it.  This  process  amounts 
to  a  transfer  of  energy  from  the  gas  to  the  solid  as  pv  work. 

The  mixture  of  solid  and  gaseous  materials  which  comprises  the 
partially  reacted  material  is  assumed  to  be  represented  by  a  collection 
of  units  similar  to  that  shown  on  figure  1,  the  walls  of  the  bomb,  in 
each  case  consisting  of  surrounding  neighbors.  These  units  move  as 
entities  in  response  to  hydrodynamic  forces,  with  the  effects  discussed 
above  superimposed.  Treatment  of  the  overall  process  is  therefore 
approached  in  the  following  way. 

Changes  in  volume  of  the  gaseous  phase,  dv  ,  are  taken  as  the 
sum  of  three  effects:  Hydrodynamics  expansions  or  extraction,  transfer 
of  material  from  the  solid  to  the  gaseous  phase  and  addition  of  heat  to 
the  gaseous  phase.  Considering  the  first  and  second  of  these  together, 
one  has: 

First , 


de  =  -  pdv 
s  r  s 


From  the  solid  equation  of  state: 


(p  +S)vs  =  (7s  -  1)  eg  +  e 


(P  +/3)dvs  +  vg  dp  =  (7s  -  1)  deg  +  egdyg 


Then: 


es  d7s  0 

-  +  - +7sj  Pdvg  -  -  vg  dp 


p  dv  p 

r  s  r 


dinp  „  . 

din  v 

s 

And  for  the  gas  (/B  »  0) 
dlnp 


din  v 


g 


e  ch 

y  m  S  — £ 

g  p  dv 


(43) 


(18) 

-(7S-1)  Pdvg 


(44) 


(45) 
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Also: 


(1  -  f  )  V  +  f  V 

s'  g  s  s 


(l-f)dv  -  vdf  +  f  dv  +vdf 
s'  g  g  s  s  s  ss 


dv 


From  (44)  and  (45) 


dv 

»  _S'_s 
dv 


.  e  d'/  , 

V  [  /  +  JB/  p  -  —  -r-M 

gl  s  H  _E-dv& I 


Vs  p 


g  P  dVo 


I 


(46) 


Combining  and  solving  for  dv^: 


dv-  - 


dv  (v  -  v  )  df 
8  8  S 


1  -  f 


1  - 


v8  |Yg_f*  £&  | 


p  dvg/ 


V  (V  +/3/p-~-p£) 
g  s  p  dvs 


Treatment  of  the  third  effect  causing  change  in  vg  involves 
the  assumption  that  all  heat  due  to  chemical  reaction  remains  in  the  gas 
phase.  One  then  has: 


f  de  *  dE  -  i  pdv 
8  g  8  g 


f  de  =  pf  dv 
s  s  r  s  s 


From  the  equations  of  state: 


Also: 


de 

s 

“7s 

de 

g 

i 

-i 

stir 

1  -fs 

GV 

s 

f 

(p  +  /3)  dvg  +  vg  dp 


dv_ 


(48) 

(49) 

(50) 

(51) 

(52) 
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Elimination  of  dp,  de  ,  de  and  dv  between  equations  48-52  gives 

s  g  s 


dv 


O'  -1)  dE 

..  .s _ : _ 


g  2 

v  (i  -  f  r 

(i  -  f  )  p  7  +  (P7S  +  0)  ^  - T2— 

s  s 


(53) 


One  then  has 


dv  =  dv  1  +  dv  " 
g  g  g 


(54) 


By  similar  methods,  the  expression  for  dvg  is  found  to  be: 


dv 


dv  *  (v  -  v  )  df 
s  r'  s 


v 

f  +  (1  -  f  )  -& 
s  s  v 

s 


e  dv 

n.  +  u p-i  *r 


7g  . 


P  dv 


(55) 


and  finally: 


de  =  -pdv 

C  r  Q 


Equations  54,  55  and  56  are  integrated  through  the  problem  and  provide 
knowledge  of  eg  at  every  point  for  use  in  equations  41  and  42. 


c.  Burning  Phase  -  Composite  Materials 

The  necessity  of  a  capability  for  treating  composite  or  multi¬ 
phase  detonable  materials  requites  some  extension  of  the  grain-burning 
concept  of  reaction.  The  material  is  conceived  as  an  assemblage  of  units, 
each  consisting  of  a  spherical  grain  of  one  phase,  normally  the  oxidizer, 
surrounded  by  the  second  phase,  i.e.  binder.  The  ignition  region  is  taken 
to  occur  at  the  interface  between  the  two. 

If  both  phases  are  self-reactive,  then  the  reaction  of  the 
spherical  grain  is  described  by  an  expression  of  the  form  of  (38)  and 
reaction  of  the  binder  phase  by  an  equation  such  as  (39)  modified  to  pro¬ 
vide  for  initiation  at  some  finite  radius.  This  capability  is  available 
in  the  present  numerical  program.  However,  in  all  problems  of  this  type 
which  were  treated,  it  was  deemed  satisfactory  to  describe  the  reaction 
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of  both  phases  by  identical  expressions,  either  38  or  39,  simply  using 
appropriately  different  values  of  the  rate  parameters.  No  distinction 
is  made  in  the  physical  properties  of  the  various  solid  phases  so  that 
a  single  solid  equation  of  state  is  used  and  a  single  value  of  eg  is 
computed  in  the  integration. 

In  multiphase  systems,  the  possibility  also  exists  that  the 
products  of  reaction  in  one  phase  can  react  with  other  phases  upon  mixing. 
It  is  presumed  that  the  mixing  takes  place  more  slowly  than  the  solid 
surfaces  recede  as  a  consequence  of  the  grain  burning  or  vaporization 
reactions,  the  diffusional  reaction  then  being  regarded  as  something 
analagous  to  a  diffusion  flame  located  approximately  at  the  original 
interface  region.  It  follows  that  the  rate  of  the  process  is  approxi¬ 
mately  that  of  the  diffusion  of  material  from  a  sphere  with  the  surface 
concentration  always  zero.  From  diffusion  theory12  the  following 
approximate  expression  is  obtained  for  the  fractional  completion  of  such 
a  diffusion  process 


£4  -  l-[o.8  exp  (-  ^  t-T>[+  0.2  Pxp  |.  «^nj(57> 

This  process  is  particularly  applicable  to  the  case  of  composite 
propellants  where  the  binder  phase  is  noneelf-reactive.  However,  it 
contributes  to  the  effective  heat  release  in  the  material  only  when  oxi¬ 
dizer  particle  sizes  are  in  the  range  of  one  micron. 

The  diffusion  coefficient  D  and  the  grain  radius  Z,  both  change 
with  hydrodynamic  expansion  in  the  region  of  the  detonation  wave  behind 
the  shock  front.  If  expansion  is  significant  through  the  reaction  zones 
of  either  the  diffusion  or  grain  burning  reactions,  representation  of 
this  change  is  necessary.  It  is  much  more  likely  to  be  required  with 
the  diffusion  reaction  since  when  the  initial  value  of  z  is  greater  than 
one  micron,  as  is  normally  the  case  the  diffusion  reaction  is  very  slow. 
These  two  quantities  are  calculated  by  means  of  the  following  expressions: 

D  ■  Dg  +  Do/p  (57a) 

Z  -  Zo  (57b) 

VO 

D  is  taken  as  the  diffusion  coefficient  in  liquids,  and  Do  that  in  gases, 
a?  ambient  pressure.  ZQ  and  V  are  the  initial  grain  radius  and  specific 
volume  of  the  explosive  materlil. 
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SECTION  3 

MATHEMATICAL  PROCEDURES 


3.1  GENERAL  PRINCIPLES 


Finite  difference  methods  for  solving  systems  of  equations  such  as 
those  discussed  in  the  previous  section  are  described  by  Richtmeyer.  3  The 
concepts  and  nomenclature  used  by  this  author  are  relatively  straightforward 
and  have  been  used  in  this  work.  Many  ways  of  differencing  are  possible, 
depending  primarily  upon  the  order  in  which  the  various  needed  quantities  are 
computed,  i.e.,  upon  the  path  through  the  computational  scheme.  Not  all  of 
these  produce  satisfactory  results.  The  most  convenient  method  of  evaluating 
a  given  procedure  is  by  trial;  while  the  methods  presented  here  appear  to  be 
satisfactory,  their  present  form  is  the  result  of  continuous  modification  and 
development  throughout  the  course  of  the  program,  and  it  is  likely  that  further 
improvement  is  possible.  A  few  specific  points  merit  special  attention. 


In  the  ignition  calculation,  it  is  assumed  that  no  chemical  energy 
is  released  until  the  material  ignites,  at  which  point  all  the  energy  from  the 
ignition  process  is  added,  and  is  taken  to  be  uniformly  distributed  through 
the  charge.  This  procedure  is  based  upon  the  precept  that  during  the  initial 
stages  of  ignition,  when  reaction  is  relatively  slow,  there  will  be  time  for 
conduction  out  of  the  ignition  region.  In  the  latter  stages,  an  exponential 
increase  in  rate  occurs  so  that  most  of  the  heat  release  occurs  in  a  very 
short  period  of  time.  The  procedure  used  is  thus  a  reasonable  approximation 
to  the  behavior  which  is  to  be  expected.  It  has  the  consequence  that  ignition 
behavior  is  defined  primarily  by  hydrodynamic  heating  and  is  relatively  in¬ 
dependent  of  the  chemical  energy  release. 


The  differencing  of  Equations  49  and  50  which  appear  in  the  grain 
burning  phase,  involves  evaluation  of  the  contribution  of  chemical  reaction 
to  the  energy  density  in  the  gas  phase,  dE.  This  was  done  as  follows. 
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Burning  is  a  process  by  which  heat  and  mass  are  added  to  the  gas 
phase.  If  this  is  conceived  as  occurring  in  discrete  steps  in  which  the  mass 
is  added  first,  and  the  energy  second,  one  can  write 


dE  - 


n 

e 


g 


(58) 


However  the  overall  process  is: 


,n+l  n+1  ,  n  n  „  _ 

f  e  -  f  e  •  Q,F„ 
g  g  g  g  2 


f  n  fn+l 
2  '  2 


+  *3P3  f3  -  f3n+1 


+  V4 


f  n  n+1 
4  '  k 


+  e 


fn  _  fn+l 
s  s 


(59) 


Transposing  and  subtracting  e  nf  n+1  from  both  sides: 

g  8 


dE  = 


n+1  n  -  n+1  n  n 

e  -e  f  =  e  -e 
8  g  g  s  g 


-  CV  v2 


n  n+1 
2  "  *2 


+  Q3F3 


fn  n+1 
3  ‘  3 


+  -  C) 


(60) 


Considerable  difficulty  has  been  encountered  in  computing  the 
correct  ultimate  velocities  and  Chapman-Jouguet  characteristics  for  detonation 
waves,  as  predicted  by  steady  state  theory  and  reported  from  experimental 
observation.  Many  explanations  for  this  are  conceivable;  however,  the  one 
which  is  now  believed  to  be  correct  is  related  to  the  use  of  the  von  Neuman- 
Richtmeyer  "q."  This  term  was  originally  used  in  numerical  procedures  for 
integrating  through  shocks  in  unreactive  media.  The  basic  equations  in  such 
problems,  consist  of  the  conservation  equation  and  an  equation  of  state,  which 
together  define  shock  behavior.  It  was  shown  that  although  use  of  the  "q" 
term  causes  distortion  in  the  shock  region,  it  has  no  effect  on  quantities 
such  as  energy,  pressure,  and  particle  velocity  at  some  distance  on  either 
side  of  the  shock  region,  and  provides  correct  treatment  of  propagation 
velocity.  However,  for  detonations,  definition  of  wave  behavior  requires  an 
additional  relation  generally  known  as  the  Chapman-Jouguet  condition.  It  now 
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appears  that  under  some  circumstances  use  of  the  von  Neuman  "q"  can  lead  to 
an  Incorrect  Chapman-Jouguet  condition.  This  point  can  be  understood  by  means 
of  the  following  arguments. 

Figure  2a  shows  an  actual  pressure  profile  through  a  detonation  wave 
as  envisioned  in  terms  of  the  hydrodynamic  theory.  The  wave  is  initiated  as 
a  shock  in  which  the  pressure  rises  discontinuously  from  pQ  to  P2 •  Reactions 
start  at  pressure  P2  and  continue  until  completion  at  the  Chapman-Jouguet 
point,  pj. 


Figure  2b  shows  the  same  wave  computed  using  equations  which  include 
the  von  Neuman  "q".  Here  the  wave  initiates  with  a  gradually  rising  continuous 
pressure  front  which  has  a  rounded  off  peak  pressure,  p_  .  In  between  there 
is  a  region  of  steeply  rising  pressure.  The  "q"  used  in  the  calculation  is 
computed  continuously  through  the  rising  pressure  region  as  proportional  to 
the  square  of  the  rate  of  compression,  and  provides  for  the  same  total  energy 
deposition  at  the  peak  of  the  wave,  p^,  as  in  a  shock  process.  In  both  cases 
the  propagation  velocity  is  given  by 


U  = 


(61) 


With  the  situation  represented  as  in  Figure  2b,  it  is  evident  that 
if  chemical  reaction  in  the  charge  is  readily  initiated,  such  initiation  may 
take  place  during  the  initial  part  of  the  pressure  rise  and  reaction  may  be 
largely  complete  by  the  time  the  material  has  passed  through  the  front  and 
reached  the  point,  p£.  Thus,  in  contrast  to  the  situation  represented  in 
Figure  2a,  a  major  part  of  the  chemical  reaction  and  the  corresponding  heat 
release  occurs,  in  effect,  as  part  of  the  shock  process.  This  has  been 
observed  to  occur  in  typical  examples  of  computed  waves. 

The  consequences  of  such  behavior  in  terms  of  wave  velocity  can  be 
appreciated  by  consideration  of  a  typical  Hugoniot  diagram  as  shown  on  Figure 
3.  According  to  the  conventional  interpretation,  ’  Figure  2a  is  repre¬ 
sented  on  such  a  diagram  as  a  process  in  which  an  increment  of  material  starts 
at  0  as  undisturbed  material,  junps  to  2  on  passing  through  the  lead  shock 
and  then  moves  down  the  line  21  while  undergoing  reaction,  to  1,  the  Chapman- 
Jouguet  point.  At  point  2,  energy  deposited  by  the  shock  in  unreacted 
material  is  given  by  area  042.  In  the  subsequent  expansion  to  point  1,  work 
represented  by  area  4216  is  performed  and  the  reaction  energy,  Q,  is  released. 
Thus,  at  point  1,  energy  remaining  in  the  system  is  given  by  area  061  +  Q. 
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FIGURE  2.  WAVE  PROFILES  WITH  AND  WITHOUT  VON  NEUMAN  0 
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FIGURE  3.  HUGONIOT  DIAGRAM  FOR  DETONATION  PROCESS 
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In  the  representation  of  the  process  shown  in  Figure  2b  on  a 
Hugoniot  diagram,  it  is  convenient  to  consider  the  extreme  case  in  which  re¬ 
action  has  gone  to  completion  by  the  time  point  2  is  reached.  In  this  cir¬ 
cumstance,  points  1  and  2  coincide  and  on  Figure  2  the  system  goes  directly 
from  point  0  to  point  1.  Since  this  is  treated  as  a  shock  transition,  energy 
equivalent  to  area  061  is  deposited  in  the  process.  The  system  energy  at 
point  1,  therefore,  is  equivalent  to  area  061  +  Q,  exactly  as  before. 

The  first  process  can  be  considered  as  an  isentropic  compression 
along  the  adiabat  053,  followed  by  heating  (via  shock  dissipation)  along  32 
with  a  corresponding  increase  in  entropy  given  by  /  TdS  =  area  032,  thence 

3 

isentropic  expansion  along  the  adiabat  27,  followed  by  heating  due  to^reaction 
along  71,  the  latter  involving  a  further  entropy  increase  given  by  /  TdS  =  Q 

7 

(the  heat  of  reaction).  The  second  process  is  made  up  of  compression  along  the 
adiabat  05,  and  heating  along  51  (due  to  chemical  reaction  and  shock  dissipation) 
with  entropy  increase  corresponding  to  f ^  TdS  =  area  051.  Considering  that 

5 

temperatures  are  highest  at  the  Chapman-Jouguet  point,  it  is  apparent  from 
these  relations  that  the  entropy  increase  in  the  second  process  is  smaller  than 
the  first.  Since  the  final  energy  content  is  the  same  in  both  cases,  it  must 
be  concluded  that  the  second  process  results  in  a  higher  pressure  at  the 
Chapman-Jouguet  point  than  does  the  first.  Therefore,  the  second  process 
cannot  terminate  at  point  1,  but  must  instead  arrive  at  some  point  at  a  higher 
pressure,  such  as  8.  From  equation  (1),  the  relative  propagation  velocities 
in  the  two  cases  will  be  given  by  ratio  of  the  square  root  of  the  slopes  of  the 
lines  01  and  08.  It  thus  becomes  evident  that  the  situation  represented  by 
Figure  2b  will  correspond  to  a  higher  velocity  than  the  value  deduced  from  the 
steady  state  treatment. 

In  present  procedures,  the  foregoing  effect  is  avoided  by  arbitrarily 
holding  up  all  reactions  until  the  zone  in  question  has  passed  completely 
through  the  shock,  as  determined  by  the  condition  q  =  0.  This  is  believed  to 
be  in  accordance  with  the  hydrodynamic  concept  of  a  shock  as  a  discontinuity, 
which  implies  an  instantaneous  transition  from  ambient  pressure  to  peak  shock 
pressure,  so  that  reaction  is  initiated  at  the  point  of  peak  pressure. 

For  reasons  similar  to  the  above,  it  is  also  probable  that  small 
consistant  errors  in  computing  pressure  in  the  region  of  the  shock  will  tend 
to  produce  erroneous  ultimate  steady  state  behavior.  For  example,  a  consistent 
overshoot  in  peak  pressure  amounts  to  an  effective  addition  of  energy  to  the 
wave.  Avoidance  of  such  effects  requires  accurate  treatment  in  the  shock  zone 
which  in  turn  implies  a  requirement  of  very  fine  zoning  in  this  region.  To 
meet  this  requirement  without  having  to  accept  excessive  computing  times,  a 
procedure  was  developed  to  impose  an  ultra-fine  space  mesh  onto  the  main  space 
mesh . 
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The  fine  zone  routine  is  desi^  led  to  automatically  insert  a  fine 
zone  mesh  in  front  of  the  oncoming  shock  and  to  merge  the  fine  zones  back 
into  their  original  zones  whenever  the  reaction  rates  fall  below  certain 
specified  values. 


At  the  beginning  of  every  problem,  two  counters  are  set  which  con¬ 
trol  the  insertion  and  deletions  of  the  fine  zone  mesh.  These  two  counters, 
jb,  which  define  the  lower  boundary  of  the  fine  zone  mesh,  and  jc,  which 
defines  the  upper  boundary  of  the  fine  zone  mesh,  are  initialized  to  j  , 
where  jmax  is  the  total  number  of  points  in  the  problem  at  the  beginning. 

The  insertion  of  the  first  fine  zone  occurs  at  the  end  of  the  first  microcycle, 
before  point  jtnalc  is  accelerated.  At  this  time,  a  A  x  is  subtracted  from 
XJmax  to  define  t*le  boundary  of  the  new  fine  zone  (xf)  .  A  check  is  made  to 

insure  that  xf  is  >  +  Ax/3.  When  xf  passes  this  test,  jmax  and  jc 

are  increased  by  one  and  the  mesh  is  moved  one  slot  upward  in  memory  from 

jb  to  jmax.  Xf  is  then  inserted  in  the  slot  formerly  occupied  by  xi  and 

Jb 

indeed  is  the  new  Xj  .  The  mass  associated  with  the  previous  Xj^  is  dis¬ 
tributed  between  the  old  and  the  new  zone.  The  time  associated  with  the  new 
point  is  the  minimum  of  the  time  associated  with  the  point  x.  and  the  old 

Jb-1 

Xj  .  The  program  then  returns  to  the  main  loop  and  does  point  jb  to 

The  general  scheme  for  inserting  a  fine  zone  is  much  the  same  as 
for  the  first  one,  except  for  determining  where  and  when  it  should  be  inserted. 
The  criteria  for  the  insertion  of  a  fine  zone,  in  the  general  case,  are: 


(1) 

(2) 

(3) 


U 


r 


j 

j+i 

j 


=  o 
<  o 

= 


(or) 


J+l 


=  h 


One  other  significant  factor  to  be  accounted  for  is  the  case  when 
Xf  <  XJb  1  +  ^hen  this  occurs,  no  fine  zone  is  inserted;  one  is 

subtracted  from  jb  and  the  problem  is  allowed  to  continue. 


-29- 


AERONUTRONIC  DIVISION 


At  the  end  of  every  macrocycle,  the  fine  zone  mesh  is  checked  to 
evaluate  whether  any  of  the  fine  zones  pass  the  qualifications  to  be  re¬ 
merged  into  their  original  zones.  The  qualification  for  a  fine  zone  to  be 
re-merged  requires  that  the  reaction  rate  of  f^  must  fall  below  a  preset 
value.  In  addition,  it  is  required  that  re  mergence  must  be  continuous  and 
conti gious  from  the  ignition  end. 

An  additional  procedure  for  conserving  computing  time  was  developed 
in  the  form  of  a  so-called  "variable  time  step"  (Aeronutronic  Hop  Program). 

In  effect,  it  involves  calculation  of  a  separate  time  interval  for  each  space 
zone,  based  upon  the  rate  at  which  the  variables  are  changing  in  the  zone. 
Calculations  of  the  variables  are  carried  out  only  in  those  regions  for  which 
this  t  is  less  than  the  macrocycle  time  Atc,  an  arbitrary  time  interval  by 
which  the  main  time  frame  is  advanced.  When  all  zones  arrive  at  the  time  of 
the  main  frame,  the  latter  is  advanced  by  A tc ,  new  variables  are  computed 
for  the  inactive  points,  and  the  process  is  repeated. 

Application  of  the  Hop  and  fine  zoning  techniques  to  the  treatment 
of  detonation  behavior  is  believed  to  represent  the  best  combination  of 
accuracy  and  economy  of  computing  time  that  is  presently  available. 


3.2  THE  FLOW  SCHEME  -  ONE  DIMENSIONAL  PROGRAM 

The  difference  equations  used  in  the  one  dimensional  case  are  given 
here  in  the  form  of  a  flow  scheme  of  the  computer  program.  The  fine  zoning 
routine  is  not  included.  All  variables  are  referred  to  the  space  coordinate 
j-1/2  unless  otherwise  indicated. 


3.2.1  HYDRODYNAMIC  EQUATIONS  INCLUDING  DISSIPATIVE  q  FACTOR 


At j+1/2  =  min  (At^2  , 

A  tn+1(2  I 

j+1/2  J 

(62) 

tn  +  Atn+1/2  :  tn+1 

j  J  c 

(63) 

(a) 

if  <  , 

go  to  64 

(b) 

if  >, 

itn+I/2 
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n+1/2  Atn+1/2  +  Atn'1/2 


(  P^x^j+i/2  +  ^  ^x^j-l/2 


n  n-1/2  n  n-1/2 

j -1/2  +  q j —1/2  "  PJ+l/2  '  qj+l/2 

,  .  n-1/2 


n  n 

Pj-l/2  *j 


n+1  n+1 

p  Xj 


where  T)  ■  1,  2,  3  for  slab,  cylindrical,  or  spherical  geometry 


n+1  n  n+1/2  .  n+1/2 

Xj  *  Xj  +  Uj  At 


In  equation  64  the  q  term  is  the  von  Neuman-Richtmeyer  artificial 
viscosity  and  is  given  by: 


n-1/2 
q j  —  1/2 


n  Q  2 

P  6 


n-1/2  n-1/2 

Uj  '  Uj'l 


n  n-1/2  n-1/2 

•  *cj  “j  ‘  “j-i 


n-1/2  n-1/2 


AS  -  X  /  t-  \  n 

“i-i  0 


_n- 1 /2  =  n 
q j -1/2  ° 


The  mass  behind  the  j  coordinate  is  defined  by: 


"j-1/2 


a  j -i/2  r« 
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n 

v  j  -1/2  ' 

1/p j-1/2 

„n+l 

mi-l/2 

P j-1/2 

/.  n+lf'  _  / 

I  n+1\ 


r .  :  0 

J 


(68) 


(a)  If  <  ,  go  to  ignition  routine 

(b)  If  >  ,  go  to  burn  routine 
3.2.2  IGNITION  PHASE 


7n+1  -  7  + 

s  o 


b  , 

CM 

O 

i 

I-* 

__  n+1  + 

V 

v  ' 

/ 

n+1 

f  n  /3  n+1/2 \ 

) 

p 

P  +“V  +  q 

i  nn  / 

7  / 

'  s 

»  P  / 

v(n+l 


n  n 

\  - 1 

7n+1  -  1 
s 


3n  +  /3  —  e/y 


n  i 

n  |  r  ,  n  7s 


n+1/2  / 


rri-1 


-  q 


n+1/2 


(69) 
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n+1 

n+1  ,  o 

1  _  n+1  1  / 

^,n+l  , 

e  * 

p  +  P 

|  v  -  e  / 

P's  •  1 

(70) 


When  a  shock  is  passing  through  an  unignited  zone,  the  ignition 
reaction  calculations  are  by  passed  until  the  zone  begins  to  expand. 

For  contracting  material 


n+1 


n+1 


n+1 

p  +  6 


n+1 


1/2 


At 


n+1 


X  U 


n+1 

i 


-  X 


n+1 

iiL 


c"+1  -  4e2Au^1/2+  lcj+1/2 


(71) 


(72) 


Return  to  Main  Hydro 


(73) 


For  expanding  material: 


Define:  At  =  t"  -  (t,i)n  i.e,  At  is  the 
°  J  1  j 

difference  between  main  coordinate  time  and 

zone  local  time  (t^)n  before  any  hydro  calculations 

have  been  mf.de  j 


Set: 

y  “ 


^  At  =  At*, 

n  ,  n+1 

e  +  e 


2  -  Q, 


1  -  f. 


G  =  B,  c 


-y 


(74) 


(75) 

(76) 


G  :  o 


(77) 
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(a)  if  -  ,  f  J  —  f"+1  ,  go  to  71 


(b)  if  4  ,  continue 


Set  St 


Jl 


fjc 


(78) 


6t  i  At 


(79) 


(a)  if  >  ,  go  to  80 


(b)  if  <  ,  go  to  82 


,n+l 


,n  -G  6  t 
fl  6 


(80) 


n+1 
1 j -1/2 


:  r 


(a)  If  >  ;  At  -  6 1  -*>  At 

go  to  75 


(b)  If  <  ,  1  -  tx 

go  to  84 


n+1 


(81) 


,n+l  -G  At 

fl  “  6 


(82) 


f^1  :  r 

j-1/2 


(83) 
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(a)  If  >  ,  go  to  71 


(b)  If  <  ,  1-  f“ 

continue 


n+1 


j  J 


t.  +  At 


At 


(83) 


go  to  71 


3.2.3  BURNING  PHASE 

The  variables  are  initialized  as  follows  upon  first  entry 
into  the  Burn  Routine. 


Grain  Surface  Burning 


n  n 

e  =  e 
s 


(85) 

(86) 


e  =  e 


(87) 


1.0 


(88) 


fU3)„ 


1.0 


(89) 


=1/3  n 


1.0 


(90) 
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n  n 

V8  "  V  J  - 1/2 


(91) 


n  n 

Vg  "  Vj-l/2 


'Ay 


g/Avo 


(92) 

(93) 


Ay 


s/Av 


fl  “  F2+F3+Fr 


(94) 


(95) 


4/3 


n+1 


1  -  F, 


1/3 


(96) 


r 


!  t 


4/3 


n+1 


1  -  F, 


1/3 


Internal  Grain  Burning 

Same  as  above  with  the  following  exceptions: 


(97) 


4/3 


(98) 


4/3 


(99) 


4/3 


-  F 


1/3 
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(101) 


(102) 


i: 

i: 

r 

E 


(  ? 

; 

t 


The  burning  phase  equations  are  the  following: 


zn  -  Z° 


/  *  t'3 


\  v 


2 

7T 

D  +  D 

.  s  ( 

,/ 

n‘ 
P  g 

ft"  -  X  1 

h _ Ll 

R] 

n 

f^-1  -  0.8  exp  (-D1*1)  +  0.2  exp  [-4D1 


n+l\ 


f1*1  .  F 

r4  1 


(103) 


(104) 


(105) 


„n+l 


(a)  If  £  ;  Fx  -*  f4  ;  go  to  107 

(b)  If  >  ;  go  to  107 


—  n 
e 


I 


:  0 
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(a)  If  <  ;  fj73 


kT1 


go  to  110  or  111 


(b)  If  £  ;  go  to  108 


1/3  n  aLAt.  Pn 


m 

-2-+1 

w 


CV  A2 

R  e"  +  ^2 

yn 

8 


ffn!^ 

zn 


c  A, 
v  3 


en  +  Q3 


For  Surface  Burning: 


-  f2  ff 1  +  f3  ff1  +  fr 


For  Internal  Burning: 


f^1  -  F_  [l-f,  n+1  +  F,  l-f-  0+1  +  F  (111) 
s  2  2  J  J  r 


38 


I 

I 


(120) 


(121) 
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yr1 


7o  + 


n+1  n+1 


-1 


(124) 
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(125) 
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(126) 


nfl 


n+1  n+1 
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7 0+1  -1 
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(127) 


n+1 


en  + 
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n  n+1 \  /  n  n+1  „  n+1/2 

vs  ■  v8  HP  +  P  +  2q 


(128) 


n+1  _ 1 


n+1 


fH+1  n+1 

s  'a 


1  -  f 


n+1 


7n+l  pn+l  +  p  fn+l 


(129) 


At 


n+1/2 


x  A  x 


nfl 

A— 


n+1  ,  Q  2  .  n+1/2  ,  ,  .  n+1 

c  -  49  +  2  \  c 


(130) 
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Return  to  main  hydro 

Definitions: 


(131) 


(132) 


(133) 
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SECTION  4 

INPUT  DATA 


4.1  EQUATIONS  OF  STATE  -  GAS  PHASE 

The  only  parameter  appearing  in  the  gas  equation  of  state, 
equation  10,  is  the  7.  This  is  taken  to  be  a  function  of  volume,  defined 
by  equation  14  of  Section  2.  As  discussed  in  Section  2,  the  term  b 
appearing  in  14  is  determined  to  have  a  value  of  1.03  by  fitting  to  ex¬ 
perimental  data  obtained  by  M.  A.  Cook.  The  term  7 0 refers  to  the  specific 
heat  ratio  of  a  normal  gas  at  ambient  conditions  and  was  assigned  the  value 
of  1.3 


4.2  EQUATION  OF  STATE  -  SOLID  PHASE 

The  required  input  for  the  equation  of  state  of  the  solid  phase, 
equation  18,  are  values  or  functions  which  define  the  parameters  7 ,  £  and 
£.  An  exper  imental  program  for  procurement  of  the  necessary  information 
is  discussed  in  a  later  section.  It  has  only  recently  been  completed, 
however,  and  the  data  which  it  has  provided  have  not  become  available  in 
time  to  be  used  in  the  theoretical  calculations  on  detonation  behavior 
which  are  reported  herein.  For  these  calculations,  preliminary  estimates 
of  the  above  parameters  have  been  made  in  the  following  way. 


As  discussed  in  a  previous  section,  a  reasonable  approximation 
is  that  at  detonation  pressures,  the  equations  of  state  of  solids  and 
gases  are  similar.  It  was  accordingly  assumed  that  the  y  for  the  solid 
could  be  calculated  by  means  of  an  expression  identical  to  that  for  the 
gaseous  phase,  equation  14,  i.e., 


7 


8 


<vx> 


(14S) 
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The  term  y  in  this  equation  le  assigned  a  value  appropriate  to  a  harmonic 
oscillator°solid,  namely  1.3.  b  is  determined  from  the  data  of  M.  A.  Cook 
to  be  1.03  as  discussed  in  Section  2. 

The  quantity  /3  was  chosen  to  be  consistent  with  the  approximate 
isothermal  compressibility  shown  by  rubber- like  materials  at  ambient 
pressure.  This  was  regarded  as  reasonable  since  the  compressibility  of 
such  materials  decreases  rapidly  at  Increasing  pressure,  and  a  major  part 
of  the  energy  deposition  during  compression  will  occur  at  low  pressures. 
Data  given  by  Brldgeman,^  and  by  Weir,  for  the  compressibility  of 
rubbers  from  0  to  1000  atmospheres  are  shown  in  Table  2.  It  indicates 
an  average  value  of  .1  dv  in  the  range  of  3x10*^.  From  Equation  18  one 


JL  dv  ,  1 

v  dp  p  +|3 

so  that  at  zero  pressure 

/3  -  3  x  104  bar 

It  is  normally  to  be  expected  that  propellants  will  be  somewhat  harder 
than  rubbers,  so  that  the  value  actually  used  in  the  calculations  was: 

£  -  5  x  104  bar. 

The  parameter  e  appears  in  expression  18  primarily  to  permit  a 
finite  volume  and  pressure  at  e  -  0.  Its  value  is  calculated  from  the 
material  density  at  ambient  conditions  for  any  specific  problem. 


4.3  THE  REACTION  RATE  EQUATIONS 

The  parameters  appearing  in  the  rate  equations  are  of  two  types: 
those  which  differ  only  to  a  small  extent  among  different  materials,  and 
those  which  differ  among  different  materials  to  such  a  degree  as  to  pro¬ 
duce  significant  differences  in  rate.  Included  in  the  first  group  are 
a,  Z,  F,  w,  Dq,  Dg.  The  second  includes  m,  Ap  .  Ap  Bp  Bg,  Bp  Qp 

Q2»  $4’ 

Values  were  deduced  for  members  of  the  first  group  in  the 
following  way.  The  term,  a,  was  taken  as  the  effective  thickness  of  a 
molecular  layer,  and  assigned  the  average  value  of  4  x  10"®cm.  Dq  and 
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TABLE  2 

.  COMPRESSIBILITY  OF  RUBBERS 

Average  Compressibility  (bar  1) 


Material 

(0-1000  bar) 

17 

3.35 

-5 

Hycar  Ml- 25 

x  10 

Hycar  OR- 15 17 

3.35 

x  10'5 

„  17 

Neoprene 

2.8 

x  10'5 

Thiokol  ST17 

2.73 

x  10'5 

Gum  Rubber16 

2.6 

x  10’ 5 
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D  were  assigned  average  values  observed  for  the  diffusion  coefficients 

8  2  /  -4  2/ 

of  gaseous  and  liquid  materials,  i.e.  0.1  cm/  and  10  cm/  respect- 

/sec  /  sec 

ively.  These  values  were  used  for  all  materials  and  all  problems. 

The  terms  Z  and  F^,  are  independent  of  the  chemical  nature  of 
the  materials  but  depend  upon  its  physical  structure.  Values  of  Z  used 
in  different  problems  varied  from  3xl0"^cm,  where  diffusion  controlled 
processes  take  place  at  significant  rates,  to  0.14  cm,  the  largest  values 
of  grain  size  reported  for  TNT.  F,  was  more  or  less  arbitrarily 

assigned  values  in  the  range  of  0.001  to  0.10.  Some  information  on  pro¬ 
pellants  was  obtained  indicating  that  porosities,  to  which  F,  can  be 
directly  related  as  discussed  previously,  are  generally  in  the  range  of  1%. 

The  parameter  w  which  defines  the  way  in  which  the  grain  burning 
reaction  becomes  insensitive  to  pressure  was  arbitrarily  chosen  to  be  20 
kllobara.  Use  of  this  number  means  that  the  grain  burning  reactions  are 
only  slightly  affected  by  pressure  above  100  kilobars. 

Parameters  of  the  second  group  relate  to  the  chemical  nature 
of  the  material.  The  order  of  the  reaction,  m,  was  taken  as  1  in  all 
cases.  Heats  of  reaction  as  obtained  from  various  literature  sources 
were  used.  The  six  rate  parameters  A^,  A^,  A^,  B^,  B^,  B^,  are  the  most 

crucial  in  terms  of  the  prediction  of  detonation  behavior.  They  are  also 
present  the  greatest  difficulty  in  terms  of  evaluation  by  independent 
means  since  the  conditions  under  which  they  are  required  to  define  reac¬ 
tion  rate  are  attainable  only  by  detonation  or  shock  processes.  In  this 
program  it  has  been  necessary  to  make  use  of  some  rate  data  obtained 
under  more  accessible  conditions. 

Four  types  of  experimental  information  have  been  utilized. 

These  are: 


1.  Thermal  decomposition  data 

2.  Minimum  initiation  pressure  data 

3.  Linear  pyrolysis  rate 

4.  Differential  thermal  analysis 

IQjB 

The  first  type  refers  to  data  reported  by  M.  A.  Cook  and  used 
by  him  to  calculate  reaction  zone  widths  for  detonation  waves  in  TNT  on 
the  basis  of  the  grain  burning  concept.  The  second  refers  to  a  method 
by  which  the  rate  constants  for  the  ignition  reaction  are  adjusted  until 
computed  values  for  the  minimum  pressure  for  initiation  are  in  agreement 
with  results  reported  from  the  experimental  measursaents  of  such  pressure. 
Experimental  data  used  for  comparison  were  reported  by  Jaffee  and  others 
at  the  Naval  Ordnance  Laboratory. The  third  type  refers  to  data  on 
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ammonium  perchlorate  obtained  by  Charken,  Anderson  and  others  at  the 
Aerojet-General  Corporation  from  their  hot  plate  experiments.  *  The 
fourth  type  of  data,  obtained  by  differential  thermal  analysis  on  double 
base  propellants  was  obtained  from  the  Allegany  Ballistics  Laboratory. 

A  summary  of  the  complete  input  used  in  calculations  on  various 
types  of  explosive  or  propellant  material  is  given  in  Tables  4-7  inclu¬ 
sive.  Examples  of  computed  detonation  behavior  carried  out  with  this 
input,  and  their  comparison  with  experimental  behavior  are  discussed  in 
the  following  section. 
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SECTION  5 

CALCULATIONS  AND  DISCUSSION 


5 . 1  INTRODUCTION 

The  mathematical  techniques  discussed  in  the  previous  sections 
permit  an  examination  of  detonation  phenomena  in  a  detail  not  possible 
by  experimental  methods,  through  the  medium  of  what  may  be  called  computed 
detonation  waves.  These  are,  in  effect,  mathematical  representations  of 
actual  detonation  processes,  and  in  the  course  of  their  creation  data  on 
pressure,  specific  volume  of  reacted  and  unreacted  material,  energy 
density  of  reacted  and  unreacted  material,  particle  velocity  and  propa¬ 
gation  velocity  are  computed  at  all  points  in  the  wave  as  a  function  of 
t ime  and  propagation  distance.  Such  information  cannot  presently  be 
obtained  from  the  observation  of  experimental  detonations.  Conversely, 
however,  the  information  is  useful  only  if  there  is  confidence  that  it 
truly  represents  reality  and  this  can  arise  only  from  extensive  demon¬ 
stration  of  conformity  between  the  characteristics  of  the  computed  and 
experimental  waves.  It  is  axiomatic  that  such  demonstrations  will  have 
to  be  in  terms  of  those  characteristics  of  experimental  waves  which  are 
observable. 

The  discussion  in  this  section  will  be  concerned  with  the  gen¬ 
eral  characteristics  of  computed  waves.  It  will  not  be  attempted  to 
present  a  detailed  comparison  between  experimental  and  computed  waves 
because,  as  discussed  under  Section  4,  a  great  deal  of  uncertainty  exists 
as  to  the  input  data  which  should  be  used  in  the  calculations  for  specific 
materials.  In  fact  when  calculations  are  stated  in  the  following  "to  be 
carried  out  for  a  specific  material"  this  should  be  understood  to  mean 
only  that  the  input  data  were  selected  to  simulate  that  material  as 
closely  as  possible.  Nevertheless  pertinent  examples  of  conformity 
between  computed  and  experimental  behavior  will  be  indicated. 


-48- 


AERONUTRONIC  DIVISION 


Examples  of  Che  latter  obtained  during  the  early  work  on  the 
program  were  presented  in  a  recent  publication^  and  will  not  be  repeated 
here.  It  is  to  be  noted  that  the  mathematical  techniques  and  computer 
routines  which  are  used  have  been  subjected  to  continuous  modification 
and  improvement  throughout  the  course  of  this  program.  Consequently, 
the  early  calculations  may  differ  in  detail  from  some  of  those  performed 
more  recently.  However,  in  no  case  do  the  later  results  requires  signifi 
cant  alteration  in  the  conclusions  derived  from  the  earlier  work. 


5.2  GENERAL  CHARACTERISTICS  OF  COMPUTED  WAVES 

The  calculations  provide  wave  trajectories,  a  typical  example  of 
which  is  shown  in  Figure  4  Input  for  this  case  corresponded  to  TNT  and 
was  as  given  in  Table  4  with  the  difference  that  Z  »  0.035  It  is  noted 
that  the  final  steady  state  velocity  attained  was  11.8  mm/i^  sec  which  sub¬ 
stantially  exceeds  the  theoretical  value.  This  problem  was  discussed  in 
Section  3  and  will  be  considered  further  in  the  following. 

Typical  curves  showing  the  reaction  and  pressure  profiles 
through  a  wave  in  TNT  with  Z  ■  0  14  cm  are  6hown  in  Figure  5,  and  profiles 
of  v  ,  v  and  e  are  shown  in  Figure  6.  Input  data  for  these  latter  two 

6  8  S 

curves  was  as  shewn  in  Table  4 

Similar  calculations  have  been  carried  out  on  composite  materials 
and  composite  double  base  materials  with  input  data  taken  as  shown  in 
Tables  3  and  5 

Figure  7  gives  the  velocity  of  a  wave  in  a  polyurethane  -APC 
propellant  initiated  with  a  100  kbar  shock.  An  initial  velocity  of  4.2 
mm/jusec  is  obtained  which  subsequently  decays  and  finally  the  wave  fades 
completely,  after  progressing  a  distance  of  about  10  centimeters  into  the 
charge.  The  explanation  for  this  behavior  is  shown  by  the  reaction-zone 
profile  of  Figure  8,  in  which  the  degree  of  completion  of  the  three  main 
reaction  processes,  f^,  f f^,  is  plotted  as  a  function  of  position  in 

the  charge  and  at  a  time  of  19  4asec  after  Initiation  It  is  seen  at 
once  that  at  this  point,  only  the  ignition  reaction,  f^,  is  proceeding 

to  any  appreciable  extent.  The  energy  derived  therefrom  is  evidently 
not  sufficient  to  maintain  the  wave,  and  decay  ultimately  ensues 

Similare  profiles  through  a  wave  in  a  composite  double  base 
material  are  shown  in  Figure  9,  This  figure  is  of  interest  in  that  it 
shows  only  slight  reaction  of  the  ammonium  perchlorate.  The  binder 
material,  on  the  other  hand  reacted  completely  at  first  in 
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FIGURE  4.  EFFECT  OF  PARAMETER  b  ON  WAVE  TRAJECTORY 
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TABLE  3 

INPUT  DATA  FOR  POLYURETHANE-APC  PROPELLANT 


-  50  kilobsrs 
*  466  cal/gm 
Bx  -  1.6  x  I012 

B2  *  8.8  *  1013 

B3  -  0 

Ax/R  -  5500°K 

A2  R  -  11000°K 

A3/R  «  11000°K 

Z  *  3  x  10  3cm 

m  “  1 

t  ■  3/iSec 
P  ' 


1.6 

0.10 

0.75 

0.15 

0 

1000  cal/gm 

380  cal/gm 

-76  cal/gm 

696  cal/gm 
0.02  cal/gm-deg 
100  kbar 


VR 

A../R 


”1 


TABLE  4 


795  cal/gm 
50  kbar 
1130  cal/gm 

1130  cal/gm 

0 

0 

4715°K 
9430°K 
1.08  x  1013 


INPUT  DATA  FOR  TNT 


B 

B 

F 


2 

3 

1 


F 


2 


F 


3 


F 


1 

Z 


a 

P 


c 

p 


1.08  x  in1 

o 

o.oi 

0.99 

0 

0 

0. 14  cm 
1.5  gm/cm3 
100  kbar 
1.5  Msec 
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FIGURE  5.  REACTION  AND  PRESSURE  PROFILES  THROUGH  COMPUTED  DETONATION  WAVE  IN  TNT. 


FIGURE  6.  SOLID-PHASE  ENERGY  DENSITY  AND  SOLID  AND  GAS  SPECIFIC 
VOLUMES  THROUGH  A  DETONATION  WAVE  IN  TNT 
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TABLE  5 

INPUT  DATA  FOR  COMPOSITE  DOUBLE 

t  »  664  cal/gm 
/3  *  50  kbar 
•  622  cal/gm 

f 

Q2  *  Q2  *  380  cal/gm 

I 

Q3  *  Qj  *  864  cal/gm 

Q4  -  0 

Aj/4  «  5500°K 
A2/4  -  11  000°K 
A3/4  «  19,400°K 
bl  =  1.6  x  1012 


BASE  PROPELLANT 

B  -  1.6  x  10 12 

2  15 

B3  -  4.6  x  lO13 

f3  -  0.01 

F2  *  0.205 

F3  -  0.470 

F .  «  0.315 

Zj  *  2.5  x  10  ^  cm 

Z3  »  12  x  10  3  cm 
2 

p  *  1.80  gm/cm 
a 

P  -  100  kbar 
t  »4  Msec 


e 

0 


At/R 

a2/r 

a3/r 


Q 

Q 


l 

1 

2 


TABLE  6 

INPUT  DATA  FOR  DOUBLE-BASE 

664  cal/gm 
50  kbar 

4.6  x  1015 

0 

0 

9700°K 
19400°K 
19400°K 
864  cal/gm 
Q2  *■-  864  cal/gm 


PROPELLANT 
Q;  -  Q3  -  0 

Q4  “  o 

Z1  -  2.5  x  10-2 

3 

Pa  »  1.80  gm/cm 

f3  -  0.001 

F2  -  0.999 
F3  -0 
Ft  -  0 

P  ■  100  kl lobar s 
tp  ■  3.0  fi sec 
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response  to  the  initial  transient.  This  reaction  was  not  maintained, 
however,  and  after  propagating  about  7cm  into  the  charge,  all  reaction 
except  the  ignition  reactions  had  ceased. 

The  explanation  for  this  behavior  appears  to  reside  in  the  high 
activation  energy,  (A^),  and  thus  the  high  temperature  coefficient  for 

the  binder  reaction.  Reaction  at  the  charge  surface  is  rapid  because  of 
the  high  energy  deposition  in  the  material  by  the  initial  shock.  However, 
the  energy  released  by  binder  reaction  is  not  sufficient  to  drive  a 
shock  of  this  intensity  so  that  following  release  of  the  surface  pres¬ 
sure  the  wave  decays.  As  a  consequence  of  the  decreasing  pressure  in 
the  wave,  the  energy  in  the  solid  drops  to  a  value  at  which  the  binder 
no  longer  reacts  at  a  significant  rate.  Similar  behavior  is  not  observed 
on  the  part  of  the  ammonium  perchlorate  reaction  because  its  lower  acti¬ 
vation  energy  makes  it  less  sensitive  to  changes  in  solid  energy  density. 
Following  the  point  at  which  the  grain  burning  reactions  subside,  the 
wave  is  supported  only  by  the  ignition  reaction,  and  this  constitutes 
the  observed  steady  state  process. 


5.3  EFFECT  OF  GRAIN  SURFACE  BURNING  VS  POINT  BURNING  ON  WAVE  BEHAVIOR 

Considerable  attention  has  been  given  during  this  program  to 
the  problem  of  differences  in  detonation  behavior  corresponding  to  the 
point  and  surface  ignition  models. 

Figure  10  shows  wave  trajectories  and  peak  pressures  for  a 
wave  in  material  corresponding  to  TNT,  computed  with  point  Ignition; 
otherwise  the  input  data  were  as  shown  in  Table  4.  Results  are  to  be 
compared  with  Figure  11,  computed  with  the  same  input  except  with  surface 
ignition.  Whereas  Figure  11  shews  a  propagating  wave  which  would 
undoubtedly  have  grown  to  a  steady  state,  the  calculation  represented  by 
Figure  10  produced  a  fading  wave. 

Figure  12  shows  the  reaction  profile  through  the  wave  of 
Figure  10,  and  it  is  of  interest  to  compare  this  profile  with  that  shown 
on  Figure  13  corresponding  to  the  trajectory  of  Figure  11.  Considerable 
difference  between  them  is  expected  due  to  the  fact  that  the  surface 
process  starts  at  a  maximum  rate  and  deaccelerates ,  whereas  the  point 
process  starts  at  a  minimum  rate  and  accelerates.  The  actual  difference 
is  substantial,  as  is  shown  in  the  two  figures. 

In  order  to  study  the  character  of  both  types  of  process  over 
the  entire  growth  period  from  initiation  to  steady  state,  the  calculations 
were  repeated  with  a  smaller  grain  size  to  produce  a  stable  wave  and  more 
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PROPAGATION  OF  COMPUTED  DETONATION  WAVE  IN  TNT 
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FIGURE  12.  REACTION  AND  PRESSURE  PROFILES  THROUGH  COMPUTED  DETONATION  WAVE  IN  TNT. 


FIGURE  13.  REACTION  ZONE  PROFILE  THROUGH  DETONATION  WAVE  IN  TNT 
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rapid  acceleration  in  both  cases.  Using  a  value  of  Z  ■  0.035  cm,  the 
trajectories  shown  in  Figures  14  and  15  were  obtained.  Also  shown  are 
wave  peak  pressures.  The  scatter  in  the  points  indicating  the  latter  is 
due  to  oscillation  in  the  numerical  integration  which  is  especially 
severe  in  the  region  of  peak  pressure.  Since  the  period  of  oscillation 
is  a  few  time  cycles  (  0.02  fisec)  which  is  short  relative  to  the  time 
required  for  changes  to  occur,  it  is  proper  to  represent  the  pressure  in 
this  region  by  the  appropriate  average  value. 

The  curves  show  pronounced  difference  in  the  growth  behavior 
for  the  two  ignition  models.  Point  ignition  produces  a  wave  which,  except 
for  the  period  during  which  it  is  overdriven  by  the  initial  pulse,  propa¬ 
gates  for  about  15  seconds  at  a  relatively  low  velocity  of  4.0  nsn/juisec 
before  accelerating  to  its  final  value  of  11.6  mm/jusec  at  23  /! sec.  Sur¬ 
face  ignition  leads  to  a  wave  which  accelerates  rapidly  from  initiation, 
attaining  final  velocity  in  about  8  fi sec. 

Figure  16  shows  reaction  zone  profiles  through  the  two  waves 
at  three  points  in  their  history.  These  generally  conform  to  expecta¬ 
tions.  The  decrease  in  reaction  zone  widths  with  time  is  due  to  the  in¬ 
crease  in  energy  deposition  by  the  lead  shock  with  increasing  propagation 
velocities,  which  results  in  increased  values  of  e  . 

The  foregoing  results  on  wave  growth  behavior  are  believed  to 
Indicate  that  the  point  ignition  process  produces  a  wave  conforming  more 
closely  to  experimental  behavior,  since  much  experimental  data  on  wave 
trajectories  follow  the  general  pattern  of  a  slow  initial  period  followed 
by  a  relatively  fast  transition  to  a  high  velocity,  essentially  as  shown 
by  Figure  14.  However,  the  trajectory  shown  on  Figure  11  for  a 
wave  computed  with  surface  ignition  but  with  a  larger  value  of  Z  also 
shows  a  low  velocity  regime,  which  is  seen  to  extend  for  more  than  four¬ 
teen  microseconds.  Therefore,  the  comparison  of  experimental  and  computed 
wave  trajectories  cannot  reliably  serve  to  distinguish  between  the  two 
types  of  Ignition  unless  independent  Information  as  to  the  appropriate 
value  of  Z  is  available.  It  is  worthy  of  note  here  that  if  loss  terms, 
such  as  might  result  from  lateral  expansion  in  the  case  of  finite  charge 
geometry,  were  included  in  the  formulation  of  the  problem,  it  is  entirely 
possible  that  with  either  type  of  ignition  the  low  velocity  period  could 
be  indefinitely  stabilized.  It  would  then  be  equivalent  to  the  phenomena 
commonly  known  as  low  order  detonation. 


FIGURE  15.  PROPAGATION  OF  COMPUTED  DETONATION  WAVE  IN  TNT 
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5.4  CALCULATIONS  OF  MINIMUM  PRESSURE  FOR  INITIATION 

As  discussed  in  Section  4,  a  factor  of  particular  importance 
in  calculating  the  detonation  behavior  of  real  materials  is  the  selection 
of  proper  values  for  the  rate  parameters  A^,  A2>  and  B2>  A  particu¬ 
larly  interesting  source  of  such  data  for  1NT  is  in  some  experimental 
work  carried  out  at  the  Naval  Ordnance  Laboratory  on  the  shock  pressures 
necessary  for  initiation  of  detonation  in  TUT. 0  Such  pressures  are 
strongly  dependent  upon  the  values  used  for  B^  and  A^  and  the  data  can 

therefore  be  used  to  evaluate  these  parameters  if  the  other  needed  infor¬ 
mation  about  the  material  is  available.  As  has  been  discussed  in  Section 
2,  the  assumption  that  the  ignition  zone  is  a  high- temperature  region  in 
the  vicinity  of  a  bubble  or  void  which  has  been  passed  over  by  a  shock, 
leads  to  the  postulate 

Ax  "  A2/2 


The  procedure  for  obtaining  values  of  the  rate  parameters  is 
to  compare  computed  values  of  the  minimum  pressure  for  initiation  of 
detonation  obtained  with  different  values  of  A.  and  B. ,  with  the  experi¬ 
mental  value.  The  computed  values  are  obtained  by  setting  and  B^ 

equal  to  zero  so  that  the  grain  burning  reactions  are  suppressed,  and 
by  using  a  low  value  of  F^  so  that  the  ignition  reaction  alone  is  not 
capable  of  supporting  a  detonation  wave.  The  calculation  is  initiated 
by  an  imposed  shock  of  ~  100  kbars  having  a  duration  of  1.5  jusec,  which 
is  heavy  enough  to  cause  Ignition.  The  disturbance  then  propagates  into 
the  material  as  a  fading  wave,  and  the  maximum  wave  pressure  at  the  last 
zone  which  ignites  is  the  minimum  pressure  for  initiation. 

Figure  17  shows  the  trajectory  of  such  a  fading  wave  obtained 
with  values  of  A^  and  derived  from  the  data  of  Cook19  which  leads  to 
the  values 

A2  “  43,400  cal/mol 
B2  -  1.08  x  1013 
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FIGURE  17.  POSITION  AND  MAGNITUDE  OF  PEAK  PRESSURE  IN  FADING 


whereupon 

-  21,700  cal/mol 
Bx  -  1.08  x  1013 

For  purposes  of  the  calculation  B2  and  B^  were  set  equal  to  zero.  Data 

on  the  heat  of  explosion,  density,  and  specific  heat  of  TNT  were  also 
obtained  from  Cook. 19  Values  of  the  other  Input  parameters  were  obtained 
as  discussed  in  Section  4. 

Comparison  of  the  fade  point  with  the  curve  of  wave  pressure 
shows  that  fading  occurs  at  a  pressure  of  ~90  kbars.  This  value  is  some¬ 
what  uncertain  because  of  the  rather  wide  fluctuations  which  occur  in 
the  computed  maximum  pressure.  Accuracy  could  be  increased  to  any  desired 
value  by  increasing  the  stability  conditions  imposed  on  the  integration 
and  decreasing  the  mesh  size.  This  would,  however,  increase  the  running 
time,  and  was  not  deemed  to  be  of  importance  to  these  calculations.  The 
experimental  value  reported  by  NOL  is  ataut  30  kbar. 

It  is  deduced  that  either  the  rate  parameters  obtained  from 
thermal  decomposition  data  predict  reaction  rates  much  lower  than  are 
actually  obtained  in  a  detonation  process,  or  that  the  assumed  equation 
of  state  parameters  lead  to  a  value  of  energy  density  behind  the  lead 
shock  much  lower  than  that  actually  obtained.  At  the  present  time  it  is 
not  possible  to  choose  between  these  two  explanations.  However,  if  the 
value  taken  for  A^  is  9500,  then  the  results  are  as  shown  in  Figure  18. 

The  minimum  shock  pressure  for  initiation  here  agrees  very  well  with  the 
experimental  value.  This  A^  corresponds  to  an  A2  for  TNT  of~  20  kcal/mol 

which  is  not  regarded  as  unreasonable. 

Calculation  has  also  been  carried  out  on  the  propagation  of  a 
detonation  wave  through  a  representative  composite  double-base  propellant 
using  the  input  data  indicated  in  Table  5.  In  this  problem,  two  grain 
burning  reactions  were  used,  one  for  the  ammonium  perchlorate  oxidizer, 
and  another  for  the  double-base  binder.  The  rate  parameters  for  the 
binder  were  obtained  from  thermal  decomposition  and  adiabatic  self¬ 
heating  data  provided  by  the  Hercules  Powder  Company.  Those  for  ammonium 
perchlorate  were  obtained  from  data  published  by  the  Aerojet  General 
Corporation  on  the  regression  rates  of  APC  surfaces  in  their  hot  plate 
experiments.  1  The  aluminum  was  treated  as  inert.  Grain  size  taken  for 
the  ammonium  perchlorate  was  the  average  size  of  the  material  used  in 
the  original  mix;  that  for  the  binder  was  1/2  the  average  distance  be¬ 
tween  APC  grains. 


FIGURE  18.  POSITION  AND  MAGNITUDE  OF  PEAK  PRESSURE  IN  FADING  WAVE  IN 
TNT  (A,/R  -  4715°K  -  A,/R  -  9430°  K) 
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The  rate  parameters  for  the  ignition  reaction  were  taken  as 
those  of  the  ammonium  perchlorate  for  two  reasons.  First,  under  condi¬ 
tions  of  low  shock  intensity  the  APC  rate  law  gives  a  faster  reaction 
than  that  for  the  binder.  This  rate  would  therefore  be  most  important 
in  the  ignition  process.  Second,  NOL  data  from  card  gap  experiments 
24,25,26  support  this  conclusion.  The  heat  of  the  ignition  reaction 
was  taken  as  the  average  of  the  heats  of  reaction  of  the  APC  and  binder 
since  they  are  assumed  to  contribute  equally  to  the  ignition  reaction. 

Figure  19  shows  the  trajectory  and  peak  pressure  in  a  computed 
wave  initiated  with  a  100  kbar  shock.  In  this  case  the  grain  burning 
reactions  were  not  suppressed  as  was  done  with  TNT  above.  Initial  wave 
velocity  is  found  to  be  about  5  mm//isec  slowing  to  about  4*'*jn/ju sec  after 
10/isec.  Peak  pressures  decay  rapidly  from  the  initial  imposed  value  of 
100  kbar  to  about  50  kbar.  Indications  are  that  this  wave  approached  a 
steady  state  condition.  However,  at  the  point  at  which  the  calculations 
were  terminated,  the  pressure  and  velocity  still  seemed  to  be  decreasing 
slightly,  and  it  is  possible  that  fading  would  ultimately  have  occurred. 

The  conclusion  is  that  the  material  is  marginal  for  detonability.  It  is 
of  considerable  interest  here  that  experimental  data  on  similar  propellant 
materials  obtained  at  the  Naval  Ordnance  Laboratory  are  consistent  with 
the  conclusions  to  be  drawn  from  these  calculations  that  a  pressure  of 
50  kbar  is  somewhat  above  the  maximum  pressure  required  for  lnltlationT  ’ 

An  important  difference  between  these  results  and  those  discussed 
previously  for  TNT  lies  in  the  fact  that  the  rate  data  used  for  the  igni¬ 
tion  reaction  in  the  present  case  were  obtained  by  the  hot  plate  technique, 
while  those  for  TNT  were  derived  from  thermal  decomposition  data.  This 
point  may  shed  some  light  on  the  question  of  which  experimental  methods 
are  most  suited  for  obtaining  kinetic  information  appropriate  for  predict¬ 
ing  detonation  behavior. 

Calculation  of  the  minimum  shock  pressure  for  initiation  of  a 
homogeneous  double-base  propellant  was  carried  out  using  the  same  tech¬ 
nique.  Chemical  and  physical  properties  were  taken  to  be  in  accord  with 
data  provided  by  the  Hercules  Powder  Company  and  in  this  case,  as  with 
TNT,  the  grain  burning  reactions  were  suppressed.  Since  this  material 
contains  no  ammonium  perchlorate,  the  data  for  the  binder  was  used  in  the 
ignition  rate  equation,  in  contrast  to  the  procedures  for  composite 
double-base  propellant.  Complete  input  data  for  the  calculations  are 
shown  in  Table  6. 

Results  of  the  calculation  are  shown  in  Figure  20,  as  plots  of 
wave  position  and  maximum  pressure  in  the  wave,  as  a  function  of  time. 

The  last  zone  ignited,  and  the  pressure  in  the  wave  upon  passing  over  this 
zone,  is  shown  by  the  vertical  line  A-A.  According  to  the  previous  dis¬ 
cussion,  this  is  the  fade  point  of  the  wave  which  thus  occurs  at  a  pressure 
of  85  kbar. 
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FIGURE  19.  TRAJECTORY  AND  PEAK  PRESSURES  OF  DETONATION  WAVE 
IN  COMPOSITE  DOUBLE -BASE  PROPELLANT 
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This  cslculstion  Is  comparable  with  experimental  determinations 
of  the  minimum  initiation  pressure  for  such  propellants  as  carried  out  by 
the  Naval  Ordnance  Laboratory.27  The  agreement  between  experimental  and 
computed  values  is  good,  and  lends  considerable  support  to  the  physical 
model  and  mathematical  techniques  which  have  been  used  in  this  program. 
However,  some  reservations  are  necessary  because  of  discrepancies, 
described  below,  which  have  appeared  in  the  properties  of  the  computed 
steady  state  wave. 


5.5  CALCULATION  OF  STEADY  STATE  CHARACTERISTICS 

When  the  basic  equations  of  a  detonation  process  are  expressed 
in  time  independent  form,  they  are  amenable  to  analytic  solution  through 
the  use  of  the  Chapman- Jouguet  condition.  The  solution  provides  a  theor¬ 
etical  steady  state  value  for  the  propagation  velocity  and  for  a  set  of 
values  of  the  pressure,  and  temperature  in  the  burnt  gas  at  the  so-called 
Chapman- Jouguet  point,  which  is  a  point  in  the  detonation  wave  defined 
by  the  condition  that  reaction  is  complete  and  the  particle  velocity 
relative  to  the  wave  shock  front  equals  the  sound  speed  at  that  point. 

An  important  test  of  the  validity  of  the  mathematical  procedures  devel¬ 
oped  under  this  program  for  computing  nonsteady  detonation  waves  is  in 
whether  or  not  the  computed  waves  ultimately  attain  the  characteristics 
predicted  by  the  steady  state  theory. 

This  question  has  been  discussed  to  some  extent  in  Section  3 
and  at  other  places  in  this  report.  It  has  been  pointed  out  that  the 
computed  waves  have  persistently  reached  steady  state  values  in  which 
the  propagation  velocities  are  much  higher  than  predicted  from  steady 
state  theory.  Such  behavior  is  shewn,  for  example,  in  Figures  4  (curve 
a),  14,  and  15,  in  which  final  velocities  in  the  range  of  11  mm//isec 
are  indicated.  This  value  compares  to  a  predicted  Chapman- Jouguet 
velocity  for  TNT  of  about  7  mm/|isec.  Pressures  at  the  Chapman- Jouguet 
point  in  these  waves  was  in  the  vicinity  of  300-400  kilobars  compared 
to  theoretical  value  of  about  100  kbars. 

The  most  recent  work  on  this  program  has  been  extensively 
concerned  with  attempts  to  ascertain  the  cause  of  this  discrepancy. 

IVo  possibilities  exist.  First,  it  may  be  due  to  an  error  in  the  input 
data  which  is  being  used.  Such  error  would  most  likely  be  in  the  values 
of  the  equation  of  state  parameters  used  either  for  the  solid  or  burnt 
gas  phases.  In  an  attempt  to  ascertain  the  sensitivity  of  the  steady 
state  velocity  to  these  parameters,  one  calculation  was  run  using  a 
value  of  0.66  for  b,  which  defines  both  7  and  7  according  to  Equa¬ 
tions  14  and  31.  Results  are  shown  in  the  two  carves  of  Figure  4 
which  were  obtained  with  identical  input  except  for  the  indicated 


-74- 


AERONUTRONIC  DIVISION 


difference  In  b.  They  show  that  use  of  b  «  0.66  gave  final  propagation 
velocities  much  closer  to  that  predicted  from  steady  state  theory  than 
were  obtained  with  the  previous  value  of  b  -  1.03. 

The  considerations  presented  in  Section  2  which  led  to  the 
assignment  of  b  -  1.03  are  not  sufficiently  unequivocal  to  justify  dis¬ 
carding  this  explanation  of  the  discrepancy  in  steady  state  wave  behavior. 
However,  together  with  additional  evidence  found  in  the  agreement  between 
computed  and  experimental  values  of  the  minimum  pressures  for  initiation, 
discussed  previously,  they  constitute  a  persuasive  argument  in  favor  of 
the  value,  1.03.  For  this  reason,  as  well  as  in  the  interest  of  confidence, 
a  re-examination  of  the  mathematical  procedures  was  undertaken  with  partic¬ 
ular  attention  to  the  possibility  that  small  consistent  errors  could 
accumulate  in  the  integration  of  the  basic  equations  and  result  in  the 
observed  discrepancy  in  the  steady  state  behavior.  A  consequence  was  the 
finding  that  under  some  conditions  the  von  Neuman  q  could  lead  to  error 
as  discussed  in  Section  3.  This  led  to  the  change  in  procedure  to  provide 
for  a  delay  in  the  ignition  until  the  material  had  passed  completely 
through  the  shock. 

The  effect  of  this  change  is  shewn  by  the  two  curves  of  Figure 
21.  Curve  1  is  the  trajectory  shown  on  Figure  15,  Curve  2  was  computed 
with  identical  input,  but  with  the  delayed  ignition.  The  pronounced 
difference  in  the  two  curves  is  an  indication  of  the  importance  of  this 
effect. 


An  additional  consequence  of  the  delayed  ignition  routine  was 
the  discovery  that  the  computational  procedure  was  depositing  shock  energy 
in  the  unreacted  solid  in  substantially  lower  amount  than  required  by  the 
Hugoniot  relation: 


Some  additional  alterations  have  been  made  as  a  consequence  of  this  dis¬ 
covery  primarily  in  the  ignition  routine  and  in  the  handling  of  the  hydro¬ 
dynamics  in  the  region  of  the  shock.  The  most  recent  modification  has 
been  the  insertion  of  the  fine  zone  routine  as  discussed  in  Section  3. 

It  is  believed  that  this  provides  a  substantial  improvement  in  the  accuracy 
of  the  treatment.  Hcwever,  as  of  the  present  date,  success  has  not  been 
achieved  in  providing  for  the  correct  deposition  of  shock  energy. 

A  wave  trajectory  run  with  the  computational  program  incorpor¬ 
ating  all  of  the  foregoing  modifications,  and  with  the  same  input  as 
Figure  21  (Table  4,  except  Z  ■  0.035  cm)is  shown  on  Figure  22.  It  is 
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believed  to  represent  the  most  accurate  treatment  presently  possible. 

It  was,  however,  not  continued  to  steady  state,  so  that  the  question 
of  computation  of  the  correct  Chapman- Jouguet  velocity  remains  unresolved. 
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FIGURE  21.  WAVE  TRAJECTORIES  SHOWING  EFFECT  OF  PROGRAMMED  REACTION  DELAY 
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SECTION  6 

TWO  DIMENSIONAL  TREATMENT 


A  number  of  Important  questions  presently  exist  concerning  the 
relationship  between  propellant  grain  geometry  and  detonation  behavior, 
and  the  effect  of  the  propellant  properties  upon  this  relationship. 

First,  there  Is  the  dependence  of  detonablllty  upon  overall 
laterlal  charge  dimensions.  For  the  case  of  a  cylindrical  charge,  this 
la  the  effect  of  charge  diameter.  One  aspect  of  this  dependence,  namely, 
the  relation  between  charge  diameter  and  steady  state  detonation  velocity, 
has  been  the  subject  of  much  of  the  classical  work  on  detonation  theory, 
and  the  ability  to  correctly  predict  the  observed  relation  constitutes 
an  Important  test  of  the  validity  of  the  theoretical  methods  devised  under 
this  program.  Other  important  aspects  of  this  problem  exist,  however, 
which  have  been  subjected  to  little  or  no  experimental  study.  Examples 
are  (1)  the  effect  of  point  initiation  vs.  initiation  by  a  plane  wave,  and 
(2)  the  effective  charge  diameter  In  the  case  of  a  cylindrical  charge 
initiated  at  a  point  half-way  between  the  two  ends.  In  all  of  such  cases 
the  most  important  item  of  information,  particularly  In  the  case  of  large 
propellant  charges,  where  experimental  data  are  lacking.  Is  the  critical 
diameter  or  dimension  at  which  a  stable  detonation  wave  can  no  longer  be 
supported. 


A  second  problem  relates  to  charges  which  have  cylindrical 
symmetry,  but  are  not  solid  cylinders.  Included  here  are  such  config¬ 
urations  as  cylinders  with  axial  cavities  of  cylindrical  or  other  shape. 
Some  experimental  data  have  been  reported^*  to  the  effect  that  such  cavi¬ 
ties  have  little  or  no  effect  and  that  only  the  outside  diameter  affects 
detonation  behavior,  at  least  until  some  critical  web  thickness  is 
reached.  Another  configuration  of  considerable  practical  Interest  is 
that  of  a  grouping  of  solid  propellant  charges  arranged  in  various 
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symmetrical  patterns  around  a  common  axis.  In  all  of  the  foregoing,  the 
intensity  and  form  of  the  initiation  event  required  to  cause  detonation, 
and  the  effect  of  location  or  other  geometrical  factors  thereon  is  of 
great  importance,  and  would  be  especially  amenable  to  investigation  by 
the  techniques  developed  during  past  work  on  this  program. 

In  order  to  provide  for  the  treatment  of  geometry  dependent 
problems  of  this  type,  a  two  dimensional  computer  program  has  been 
written  which  incorporates  the  procedures  described  in  the  foregoing 
sections.  This  program  is  based  upon  a  two-dimensional  hydrodynamic 
program  (with  cylindrical  symmetry),  employing  a  Lagrangean,  triangular 
mesh  and  incorporating  the  variable-time-step  features  introduced  in  HOP, 
the  one-dimensional  code.  This  program  was  prepared  for  Plcatlnny  Arse¬ 
nal  under  contract  DA-04-495-ORD-3095.  ®  The  original  code  associates 
two  (triangular)  zones  with  each  point  of  the  mesh.  The  triangular  zones 
give  somewhat  improved  accuracy  and  stability  than  the  more  typical  quad¬ 
rilateral  zones. 

The  extra  data  required  for  each  zone  in  order  to  Include  chem¬ 
ical  kinetics  severely  restricts  the  number  of  mesh  points  which  can  be 
internally  stored  in  a  32000  word  computer.  Accordingly,  the  triangular 
code,  called  ROC-VTS,  was  rewritten  using  a  quadrilateral  mesh  but  re¬ 
taining  the  variable  time  step  feature.  It  is  possible  to  store  about 
1200  mesh  points  internally.  Basic  concepts  used  in  this  code  are  the 
following. 


Introduce  Eulerian  coordinates  R,  the  distance  from  the  axis 
of  symmetry,  and  Z,  distance  measured  along  the  axis.  Mesh  points  are 
characterized  by  Lagrange  coordinates  K,  L  and  the  temporal  mesh  is 
characterized  by  the _ superscript  n.  The  momentum  equations  to  be  solved 
for  new  velocities,  R  and  Z  are: 


.  nti, 

*Sc,L 


hi  +  (p+q>K-%,J 


*Si,L  +|(p+q)K-^,Lf%  ^K-I.l'Vl+P  +  (p+q)K+%,L+% 

(134) 

<ZK,W-rZK+l,L)  +  (p+q)K+%,L-^ZK+l,L*ZK,L-l*  + 


(^S-%,L-%(ZJ,L-1-ZS-1,L> 


,L 
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in+ij  Ltl-k  _ 
k,  L  ZK,L 


J(p+q)K.l(iW.ij(RK.1>L-RK>I/fl)  + 
n  n  n  _  _  (135) 

'  *Sc+l ,L^  +  ^^^K+^.L-lf^+l.L'^C.L-P  + 

1 1»  '»n 

>  ,  .n  ,_n  _n  .l/Atl 


where 


At 

f* 


r 

/K, 


Atn+^  +  Atn’^ 
Q  K.L  a  K.L 


L  *pA)K-%,Lrf%  +  ^A)K+%,Lf%  +  ^A)K+%,L-%  +  (pA)K-*,L-%. 


(136) 


In  the  above,  p  Is  the  density  and  A  the  area  (of  revolution)  of  the 
quadrilateral  zone  In  question.  These  momentum  equations  are  felt  to  be 
a  substantial  Improvement  over  those  typically  used  for  two-dimensional 
cylindrical  Lagrange  hydrodynamic  codes. The  Interested  reader  Is 
referred  to  references  28,  29,  or  30  for  further  discussions  and  refer¬ 
ences  of  two-dimensional  hydrodynamic  computing  techniques. 


This  two  dimensional  program  Is  presently  running  and  has  been 
used  for  short  demonstrations.  However,  no  full  problems  have  been  run. 
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EXPERIMENTAL  INVESTIGATION  OF  EQUATION  OF  STATE 


The  objective  of  this  phase  of  the  program  was  the  determina¬ 
tion  of  the  quantities  0  and  7  in  Equation  18  as  functions  of  the  variables 
p,  v,  and  e.  The  work  involved  the  procurement  of  p-v-T  data  by  static 
pressure  techniques,  making  use  of  apparatus  and  experimental  methods 
described  by  Kennedy  for  work  in  the  ultra-high  pressure  region. 

The  experimental  method  involves  the  use  of  a  high  pressure 
piston-cylinder  press32  shown  schematically  in  Figure  23.  In  this  appa¬ 
ratus  a  tungsten  carbide  piston  is  driven  by  a  high  pressure  hydraulic 
system  into  a  tungsten  carbide  cylinder  supported  by  a  specially  designed 
pressure  plate  consisting  of  a  series  of  press-fitted  concentric  rings 
(see  Figure  24).  Pressure  within  the  cylinder  is  determined  from  ram  oil 
pressure  as  measured  on  a  calibrated  Heise  gauge,  and  the  known  ratio  of 
cylinder  diameter  to  ram  diameter.  Volume  within  the  cylinder  is  deter¬ 
mined  from  piston  position  as  measured  either  by  a  dial  gauge,  micrometer, 
or  a  linear  potentiometer,  the  latter  being  used  to  drive  one  leg  of  an 
x-y  recorder. 

Information  is  obtained  as  an  upward  and  downward  curve  of 
pressure  versus  volume  change,  first,  by  increasing  the  pressure,  in 
steps,  to  the  maximum  value  desired,  and  then  similarly  decreasing  it 
stepwise  back  to  zero,  with  measurements  of  the  piston  movement  at  each 
step.  In  all  such  curves,  the  upward  and  downward  legs  demonstrate  a 
definite  hysteresis  having  the  general  characteristics  shown  in  the 
linearized  diagram  of  Figure  25.  This  is  due  primarily  to  frictional 
forces  -  both  of  the  piston  on  the  walls  of  the  cylinder  and  internal 
friction  of  the  sample.  Thus,  for  example,  the  pressure  difference  be¬ 
tween  points  3  and  4  represents  the  vector  difference  of  the  frictional 
forces  for  the  up  and  down  strokes.  If  the  magnitudes  of  these  forces 
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FIGURE  25.  rYFTCAl  I'Oi.M  OF  PRESSURE  -  DISPLACEMENT  CURVE 
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are  Che  same  (at  the  same  displacement)  an  average  sample  pressure  at 
this  displacement  is  represented  by  a  point  half-way  between  3  and  4 
(i.e.  point  9).  It  follows  that  to  a  good  approximation  the  true  p-v 
curve  for  the  material  in -the  cylinder  is  given  by  the  curve  2,6.  Direct 
observation  of  the  curve  2,6  i*  possible  if  sliding  friction  is  relieved 
prior  to  each  observation  by  rotating  the  piston,  and  if  sufficiently 
long  times  are  allowed  between  observations  for  sample  Internal  friction 
to  anneal  out.  It  has  been  reported  by  Kennedy, ^3  however,  that  use  of 
the  mid  points  between  the  ascending  and  descending  curves  will  give  the 
true  p-v  curve  with  an  accuracy  which  is  within  the  other  experimental 
errors  of  the  experiment. 

It  is  necessary  to  convert  the  information  such  as  shown  in 
Figure  25  into  a  sample  compressibility,  and  ultimately  into  a  p-v  curve. 
To  do  so  require  s  the  input  of  certain  information  about  the  press  geom¬ 
etry  and  its  distortion  as  a  function  of  pressure.  The  latter  is  known 
as  the  press  constant. 

Figure  26  shows  a  detailed  diagram  of  the  press  set-up  during 
an  experiment.  Here  As  is  the  pressure  plate  shown  in  detail  in  Figure 
24  and  which  forms  the  high  pressure  cylinder,  and  Bs  is  the  piston. 

C8  Is  the  end  loading  platon  for  clamping  the  sides  of  the  pressure  plate 
between  the  cylinder  head  block,  Ds,  which  in  turn  rests  against  the 
press  head  Es,  and  the  spacer  Fs.  Gs  is  a  force  transmitting  column, 

H8  is  a  pillow  plate  and  Is  is  the  master  ram. 

The  sample  volume  is  defined  at  any  time  by  the  distance  be¬ 
tween  the  piston  head  and  the  cylinder  head,  i.e. 


The  change  in  this  distance  with  ram  motion  is  equal  to  the  ram  movement 
less  the  press  constant.  If  the  distances  X^  and  X8  are  known  at  any 

S  ~ 

reference  ram  position,  then  the  distance  X^  at  any  other  ram  position 
is  given  by 

X1  =  X10  +  X20  '  X30  'tAxS  +  cS  037) 

where  the  subscript  zero  refers  to  the  reference  position. 
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With  a  given  piston  diameter,  it  is  assumed  that  the  press 
constant  is  dependent  only  upon  cylinder  pressure  rather  than  ram  posi¬ 
tion,  i.e.  sample  volume.  This  is  true  only  if  the  radial  distortion  of 
the  cylinder  wall  is  independent  of  the  volume  or  other  factors  related 
to  the  sample.  Evidence  obtained  by  Kennedy  has  indicated  that  radial 
expansion  of  the  cylinder  is  negligible  under  the  pressures  considered 
here.  Consequently  if  the  sample  volume,  X^,  is  zero: 


X20  =  X30 

and  the  press  constant  is  readily  determined  as 

CS  «AXS  (138) 

g 

The  reference  point  for  measurement  of  the  ram  motion, Ax  ,  is 
generally  chosen  in  such  a  way  as  to  provide  the  smallest  value  of  the 
press  constant. 

The  foregoing  considerations  lead  to  the  following  procedure 
for  a  press  run.  First,  with  no  sample  in  the  cylinder  the  ram  is  set 
at  a  position  corresponding  to  zero  simple  chamber  volume.  This  is 
readily  recognized  as  a  point  at  which  hydraulic  pressure  rises  steeply 
with  ram  motion.  A  plot  of  ram  displacement  vs  pressure  beyond  this 
point  is,  according  to  Equation  138  also  a  plot  of  press  constant  vs 
pressure.  Little  or  no  hysteresis  is  found  in  the  up  and  down  legs  of 
this  curve. 

A  sample  is  then  inserted  in  the  pressure  chamber  and  a  press 
cycle  run  as  described  in  the  foregoing.  It  frequently  happens  that 
sliding  friction  does  not  disappear  completely  at  zero  pressure,  so  that 
the  curves  of  the  up  and  down  legs  show  unsymmetric  regions  at  both  ends 
as  indicated  in  Figure  25.  In  such  cases  the  points  in  these  regions 
cannot  be  used  for  centering.  They  are  normally  excluded  by  restricting 
the  centering  procedure  to  that  region  where  the  up  and  down  legs  are 
relatively  parallel. 

The  centered  curve  is  corrected  for  the  press  constant  and  for 
extrusion  of  the  sample  into  champfered  spaces  at  the  piston  and  cylinder 
heads.  Additional  corrections  relate  to  the  presence  in  the  cylinder  of 
items  other  than  the  material  of  interest.  These  include  a  thermocouple 
and  thermocouple  retainer  made  of  stainless  steel,  and  sometimes  a  wafer 
of  inert  material  between  the  piston  head  and  the  sample.  The  latter 
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was  used  in  working  with  combustible  materials  which  conceivably  could 
become  ignited  by  viscous  heating  if  they  were  forced  into  the  space 
between  piston  and  cylinder  wall.  Corrections  to  be  made  for  these  sub¬ 
stances  must  be  derived  from  data  on  their  compressibility  obtained  in 
separate  experiments. 

Pressure- volume  curves  at  elevated  temperatures  were  obtained 
by  clamping  ring  heaters  to  the  sides  of  the  pressure  plate  A8.  Maximum 
temperatures  investigated  were  about  100  C.  Sample  temperatures  were 
monitored  in  all  cases  by  means  of  a  chromel-alumel  thermocouple  inserted 
through  the  cylinder  head  block  and  into  a  hole  in  the  sample.  Tempera¬ 
ture  variation  over  a  cycle  was  not  found  to  be  more  than  1-2  C. 

The  p-v-T  data  obtained  are  presented  in  Figures  27  to  40, 
inclusive.  Curves  are  included  for  a  number  of  pure  or  homogeneous 
materials,  and  on  their  various  mixtures.  Tables  8  and  9  give  pertinent 
data  on  the  composition  of  the  salt-polyurethane  speciman  and  on  the  poly- 
urethane-APC  propellant  speciman.  Table  10  gives  p-v-T  data  on  the  pro¬ 
pellant  in  tabulated  form. 

A  number  of  the  figures  showing  data  on  mixtures  also  include 
curves  labeled  "calculated".  These  were  deduced  by  assuming  that  the 
specific  volumes  of  the  constituents  of  the  mixture  were  additive  at  all 
pressures.  Interest  in  this  point  derives  from  two  considerations.  In 
the  first  place,  some  mixtures  pertinent  to  this  program  are  combustible 
and  susceptible  to  ignition  at  elevated  temperatures,  whereas  their 
major  constituents  are  much  less  so  in  the  pure  state.  This  would  be 
true,  for  example,  of  a  polyurethane  ammonium  perchlorate  propellant, 
such  as  the  one  for  which  data  is  presented  in  Figure  39.  In  the  second 
place,  many  mixtures  of  interest  are  composed  of  a  relatively  small  number 
of  the  same  basic  constituents,  and  a  considerable  saving  in  time  would 
be  possible  if  their  p-v-T  relationships  could  be  derived  directly  from 
the  corresponding  curves  for  the  components. 

Agreement  between  the  experimental  and  calculated  curves  are 
fairly  good  in  some  cases,  and  less  so  in  others.  No  general  conclusions 
have  been  drawn  as  to  the  feasibility  of  this  technique.  However,  it  was 
used  to  obtain  the  p-v  curve  for  a  polyurethane  propellant  at  about  100  C. 
The  result  is  shown  on  Figure  40. 

The  information  presented  in  the  figures  can  be  used  to  evalu¬ 
ate  the  parameters  and  y  in  Equation  18,  Section  2,  in  a  straight¬ 
forward  way.  However,  this  had  not  been  done  as  of  the  date  of  this 
report. 
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FIGURE  27.  PRESSURE- VOLUME  CURVES  FOR  SODIUM  CHLORIDE  AT  24°C 
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FIGURE  28.  PRESSURE- VOLUME  CURVE  FOR  SODIUM  CHLORIDE  AT  100°C 
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FIGURE  29.  PRESSURE-VOLUME  CURVE  FOR  POLYURETHANE  AT  23°C 


AERONUTRONIC  DIVISION 
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FIGURE  30.  PRESSURE- VOLUME  CURVE  FOR  POLYURETHANE  AT  99.5  -  101°C 
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FIGURE  31.  PRESSURE-VOLUME  CURVES  FOR  AMMONIUM  PERCHLORATE  AT  21.5  -  22.06°C 
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FIGURE  32.  PRES SURE -VOLUME  CURVE  FOR  AMMONIUM  PERCHLORATE  AT  100  -  102°C 
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FIGURE  33.  PRESSURE-VOLUME  CURVES  FOR  SODIUM  CHLORIDE- POLYURETHANE 
MIXTURE  (40  -  60%  BY  WT)  AT  22  -  23°C 
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FIGURE  34.  PRESSURE-VOLUME  CURVES  FOR  SODIUM  CHLORIDE- POLYURETHANE 
MIXTURE  (40  -  60%  BY  WT)  AT  91  -  94°C 
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FIGURE  35.  PRESSURE-VOLUME  CURVES  FOR  SODIUM  CHLORIDE- POLYURETHANE 
MIXTURE  (70  -  307.  BY  WT)  AT  23*C 
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FIGURE  36.  PRESSURE- VOLUME  CURVES  FOR  SODIUM  CHLORIDE- POLYURETHANE 
MIXTURE  (70  -  30%  BY  WT)  AT  99  -  102 °C 
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FIGURE  37.  PRES SURE- VOLUME  CURVES  FOR  SODIUM  CHLORIDE-AMMONIUM 
PERCHLORATE  MIXTURE  (50  -  50%  BY  WT)  AT  21°C 
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FIGURE  38.  PRESSURE- VOLUME  CURVES  FOR  SODIUM  CHLORIDE-AMMONIUM 
PERCHLORATE  MIXTURE  (50  -  50%  BY  WT)  AT  99  -  101 °C 
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FIGURE  39.  PRESSURE- VOLUME  CURVES  FOR  ALUMINIZED  SOLID  PROPELLANT 
(B.  F.  GOODRICH  DESIGNATION  E107M)  AT  23°C 
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FIGURE  40.  PRESSURE-VOLUME  CURVE  FOR  ALUMINIZED  SOLID  PROPELLANT 
(B.  F.  GOODRICH  DESIGNATION  E107M)  AT  100°C 
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TABLE  7 

INPUT  COMMON  TO  ALL  PROBLEMS 


a 


4.0  x 


10 


-8 


cm 


7  -  1.3 
'o 

b  =1.03 

-4  2 

D  =  10  cm  /sec 

s 

-1  2 

D  =  10  cm  -bar/sec 

o 

w  =  20  kbar 


m  =  1 

Cv  =  0.3  cal/gm-deg 
TABLE  8 

PREPARATION  OF  SALT- POLYURETHANE  SPECIMENS 


Weight  %  of  Salt-Polyurethane 


Composition 

70-30% 

40-60% 

Polyproylene  glycol  2025 

24.980 

50.02 

Trimethylol  propane 

0.956 

1.91 

Ferric  acety lacetonate 

0.04 

0.07 

Toluene-2,  4-diisocyanate 

4.042 

8.00 

Sodium  Chloride 

70.000 

40.000 

Cure 

24  hrs  @  160°F 

24  hrs  <a  160°F 
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TABLE  9 

DATA  FOR  SOLID  PROPELLANT  (E107M)* 

Type:  high  energy  polyurethane  base,  castable  and  case-bondable , 
slow-burning,  composite  solid  propellant 

Theoretical  Performance:  I  (1000)  *  255  sec. 

sp 

Density  (77°F):  0.0608  lbs/in^ 

Composition: 


Components 

Weight  % 

1. 

polyurethane 

25.0 

2. 

ammonium  perchlorate 

57.3 

3. 

aluminum 

17.7 

*  "High  Energy  Composite  Solid  Rocket  Propellants  -  Polyurethane 
Base  and  "c"  -  Rubber  Base"  B,  F.  Goodrich  Aviation  Products, 
Division  of  B.  F.  Goodrich  Co.,  Rocket  Motors,  Rialto,  California. 
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TABLE  10 

COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  SPECIFIC  VOLUMES 
FOR  SOLID  PROPELLANT  (E107M  B.  F.  GOODRICH)  AT  AMBIENT  TEMPERATURES 


3 

Specific  Volumes,  cm  /g 


Pressure , 
k  bar 

Calculated 

Experimental 
(Curve  #43B) 

0 

.5939 

.5944 

4 

.5656 

.5781 

8 

.5509 

.5628 

12 

.5403 

.5502 

16 

.5315 

.5403 

20 

.5329 

.5316 

24 

.5166 

.5243 

28 

.5106 

.5178 

30 

.5076 

.5146 

Notes: 

1.  Specific  Volume  at  77  F  is  given  as  0.5945  cm  / g  by  the  manufacturer 
and  compares  closely  with  the  figures  shown  in  Table 

2.  The  experimental  and  calculated  curves  are  different  by  27„  in  regions 
of  maximum  difference. 
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